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gravitation  property,  then  density  is  a  term  having  a  precise  meaning  ; 
but  if  the  distribution  of  the  gravitation  property  through  the  space 
occupied  by  the  body  has  not  this  mathematical  continuity,  we  cannot 
attach  any  meaning  to  the  volume  integrals  till  we  have  first  invented 
a  suitable  new  meaning  for  the  term  '  density ' ;  and  the  inevitable 
vagueness  that  will  arise  in  the  new  definition  will  preserve  in  the  final 
results  that  slight  lack  of  precision  present  in  the  terms  of  statement 
of  the  gravitation  law,  which  might  at  first  sight  appear  to  have 
dropped  out  of  results  represented  by  such  precise  mathematical 
expressions  as  volume  integrals. 

It  is  practically  certain  that  no  substance  can  be  subdivided  without 
limit  into  small  portions  each  of  which  possesses  the  gravitation 
property.  There  must  be  a  stage  of  subdivision  beyond  which  the 
component  portions  cease  to  have  the  properties  of  larger  portions  of 
the  substance,  and  we  may  speak  of  the  smallest  portion  of  a  substance 
that  has  the  gravitation  property  as  a  'particle'  for  purposes  of  the  present 
discussion.  What  the  order  of  magnitude  of  a  particle  may  be  it  is 
difficult  to  guess,  but  the  kind  of  generalisation  from  large  bodies  to 
small  bodies  which  led  to  the  conception  of  an  element  of  mass  suggests 
the  possibility  that  the  process  of  subdivision  without  loss  of  the 
gravitation  property  might  be  continued  till  we  arrive  at  the  molecule 
of  Chemistry  or  Gas  Theory.  There  is  no  experimental  evidence  to 
prohibit,  and  possibly  some  to  justify  our  carrying  the  generalisation 
so  far  (provided  we  set  some  limit  to  the  smallness  of  the  distance  at 
which  the  attraction  of  two  particles  is  supposed  to  obey  the  law  of  the 
inverse  square),  and  the  great  simplicity  of  the  law  thus  obtained 
makes  it  an  interesting  one  to  study. 

If  a  body  has  not  got  mathematical  continuity  in  the  distribution 
of  the  gravitation  property  throughout  its  volume,  and  so  is  to  be 
regarded  as  made  up  of  particles,  we  may  speak  of  it  as  having 
'discontinuous'  structure.  And  if  it  be  supposed  that  the  particles 
may  be  as  small  as  molecules,  we  must  form  a  mental  picture  of  the 
structure  in  which  there  appears  no  trace  of  material  continuity,  the 
substance  being  represented  by  discrete  molecules  or  systems  occupying 
spaces  with  somewhat  indefinite  boundaries,  separated  by  more  or  less 
empty  regions  which  may  be  called  intermolecular  space;  (if  the 
molecules  move  it  may  be  supposed  that  their  motions  do  not  affect 
the  properties  under  consideration).  The  problem  of  finding  the 
potential  or  attraction  of  such  a  body  at  any  point,  if  formulated  on  a 
greatly  magnified  and  coarser  scale,  would  in  some  respects  resemble 
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the  problem  of  evaluating  the  potential  or  attraction  of  a  mass  of  sand 
or  other  granular  matter.  We  want  to  see  how  volume  integrals  present 
themselves  as  approximate  solutions  of  such  problems. 

It  is  unnecessary  to  dwell  here  upon  the  familiar  definitions  of  the 
intensity  of  force  at  a  point,  or  the  potential  at  a  point,  due  to  a 
gravitating  body.  But,  for  future  reference,  we  may  emphasize  the 
fact  that,  among  the  mathematical  properties  of  the  potential  at  a 
point  outside  the  body,  that  which  may  be  taken  as  the  fundamental 
physical  definition  is  the  fact  that  its  space-gradient  at  any  point  is 
vectorially  equal  to  the  intensity  of  force  there.  If  a  point  is  so 
situated  that  a  physical  definition  of  intensity  of  force  there  is 
impossible,  this  physical  definition  of  potential  breaks  down,  and  we 
are  at  liberty  to  substitute  some  convenient  purely  mathematical 
definition  for  which  it  may  be  possible  afterwards  to  find  a  physical 
interpretation. 

The  potential  of  a  body  of  discontinuous  structure  at  an  external 
point  P  is  the  sum  of  the  potentials  at  P  due  to  the  particles  that 
compose  the  body,  i.e.  ^mr~l,  where  m  is  the  mass  of  a  particle  and 
r  its  distance  from  P.  Here,  in  accordance  with  what  has  been  said 
above,  r  is  not  precisely  defined,  and  a  corresponding  lack  of  precision 
must  be  present  in  ^mr~l.  By  assuming  P  to  be  not  too  close  to  any 
particle  of  the  body  we  can  ensure  that  each  r  shall  always  be  great 
compared  with  the  linear  dimensions  of  the  corresponding  particle. 

When  we  endeavour  to  compare  the  values  of  this  expression  for 
the  potential  at  different  points,  we  recognise  that  the  sum  of  a  finite 
but  extremely  great  number  of  extremely  small  terms  is  a  most  trouble- 
some function  to  work  with,  and  so  there  naturally  suggests  itself  the 
device  of  getting  a  probably  very  approximate  equivalent  function  by 
replacing  the  sum  by  the  limit  to  which  it  would  tend  if  the  number 
of  terms  could  be  increased  indefinitely  while  each  separate  term 
decreased  correspondingly;  this  process  would  give  us  the  potential  in 
the  form  of  the  definite  integral  /pr~Wr,  where  dr  is  an  element  of 
volume  and  pdr  the  corresponding  mass. 

But,  as  has  been  suggested  already,  the  transition  from  a  sum  of 
terms  to  a  definite  integral  would  imply  the  possibility  of  endless 
subdivision  of  the  material  mass  into  elements  each  possessing  the 
gravitation  property,  whereas  it  is  practically  certain  that  matter 
cannot  be  so  endlessly  subdivided.  In  fact  the  use  of  the  definite 
integral  form  implies  a  regarding  of  matter  as  continuously  extended 
through  the  space  which  it  effectively  occupies,  and  attributes  to  the 
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density  p  at  any  point  the  value  obtained  by  passing  to  a  mathematical 
limit  in  the  usual  fashion,  that  is  to  say  the  limit  of  the  quotient  of 
mass  by  volume  for  a  region  surrounding  the  point  as  the  dimensions 
of  the  region  tend  to  zero.  The  molecular  view,  however,  requires  us 
to  cease  subdividing  matter  beyond  a  certain  stage,  and  so  prevents 
our  ever  arriving  at  the  kind  of  limit  which  is  known  as  an  integral. 

3.  Nevertheless  the  potential  at  the  point  P  of  an  assemblage  of 
discrete  particles  in  a  finite  region  may  be  equal  in  value  to  a  volume 
integral  taken  throughout  the  region  if  the  integral  be  supposed  to 
refer  to  a  hypothetical  continuous  medium  occupying  the  same  region 
and  having  a  suitably  chosen  density  at  each  point.  It  is  only 
necessary  to  choose  the  law  of  density  properly,  and  to  this  end  there 
suggests  itself  the  device  of  taking  for  each  point  A  some  sort  of 
average  density,  based  on  a  consideration  of  all  the  masses  within 
a  very  small  but  finite  region  surrounding  A.  The  dimensions  of  this 
small  region  might  be  settled  by  convention,  but  we  need  only  consider 
the  order  of  magnitude  of. these  dimensions. 

The  kind  of  smallness  that  we  want  in  this  connexion  is  what 
we  may  call  physical  smallness,  as  distinguished  from  mathematical 
smallness  to  which  there  is  no  limit.  Physical  smallness  implies 
smallness  which  appears  extreme  to  the  human  senses,  but  it  must  not 
be  a  smallness  so  extreme  as  to  necessitate  passage  from  molar  physics 
to  molecular  physics  ;  it  must  leave  us  at  liberty,  for  example,  to 
attribute  to  matter  occupying  a  physically  small  space  the  properties 
of  matter  in  bulk  if  these  should  be  different  from  the  properties  of 
isolated  molecules.  In  fact  a  physically  small  region,  though  extremely 
small,  must  still  be  large  enough  to  contain  a  very  great  number  of 
molecules.  It  is  estimated  that  a  gas,  at  normal  temperature  and 
pressure,  has  about  4  x  1019  molecules  per  cubic  centimetre  ;  thus  a  cube 
whose  edge  is  6  x  10~5  centimetres  (roughly  the  wave  length  of  sodium 
light)  contains  more  than  8,000,000  molecules  ;  if  we  regard  a  million 
as  a  large  number,  the  wave  length  of  sodium  light  is  (for  other  than 
optical  purposes)  physically  small,  and  it  is  known  that  very  much 
greater  lengths  than  this  appear  to  our  senses  extremely  small.  We 
are  therefore  in  a  position  to  speak  of  lengths  which,  though  extremely 
small,  are  very  great  compared  with  other  physically  small  lengths. 

Now  it  is  not  suggested  that  the  gravitation  property  is  a  molar 
property  of  matter,  not  possessed  by  a  single  molecule,  for  we  have 
adopted  just  the  opposite  hypothesis ;  and  so  it  might  be  thought 
justifiable  to  make  the  region  round  A,  used  for  calculating  p, 
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smaller  than  merely  physically  small.  This  point  will  be  referred  to 
again,  but  at  present  it  suffices  to  remark  that  the  p  generally  used  in 
potential  theory  is  a  continuous  function  whose  value  is  not  subject  to 
very  rapid  fluctuations  as  A  moves  from  one  position  to  another.  To 
get  such  continuity  and  smoothness  in  the  suggested  average,  and  to 
avoid  the  difficulty  that  would  arise  if  the  number  of  molecules  inter- 
sected by  the  boundary  of  the  region  so  as  to  be  neither  obviously 
inside  nor  obviously  outside  were  not  very  small  in  comparison  with 
the  total  number  inside,  we  must  take  account  of  a  large  number  of 
molecules  at  a  time,  and  so  we  have  to  assume  that  the  region 
surrounding  the  point  A  and  used  in  getting  a  value  for  p  is  only 
physically  small. 

As  regards  the  system  of  averaging,  it  is  clear  that  in  the  case 
under  discussion  we  get  the  best  agreement  between  the  potential  of 
the  actual  and  that  of  the  hypothetical  system  if  we  give  to  each 
molecule  an  importance  proportional  to  the  product  of  its  mass  and  the 
reciprocal  of  its  distance  from  P ;  but  it  would  be  unfortunate  to  be 
obliged  to  average  in  this  fashion,  as  we  should  thereby  get  a  value  of 
P  at  A  which  would  not  be  independent  of  the  position  of  P.  And  we 
should  get  quite  a  different  law  of  density  if  we  were  dealing  with  some 
other  integral,  say  an  attraction  integral,  instead  of  that  representing 
potential. 

So  long,  however,  as  P  is  at  a  distance  from  A  great  compared 
with  the  linear  dimensions  of  the  physically  small  region  used  for  the 
purpose  of  averaging,  the  values  of  r~l  for  the  molecules  in  this  region 
are  very  nearly  equal,  and  so  there  is  very  little  error  if  in  taking  the 
average  we  give  to  each  molecule  an  importance  simply  proportional  to 
its  mass.  We  thus  get  for  the  density  p  of  the  hypothetical  continuous 
medium  the  quotient  of  mass  by  volume  for  the  physically  small  region 
considered.  This  value  of  p  has  the  advantage  that  it  is  independent 
of  the  position  of  P,  and  of  the  particular  physical  quantity  whose 
integral  expression  is  being  investigated.  But  its  great  advantage, 
and  the  real  reason  why  we  adopt  it,  is  that  it  is  that  density  of  a 
substance  which  we  actually  arrive  at  by  practical  methods  of  measure- 
ment ;  for  ordinary  laboratory  measurings  and  weighings  are  applied 
to  portions  of  a  substance  which  are  far  from  the  limits  of  physical 
smallness,  and  so  give  us,  not  the  sizes  and  masses  of  individual  mole- 
cules, but  only  the  total  mass  and  the  total  space  effectively  occupied. 

4.  So  far  we  have  considered  only  the  case  in  which  the  point  P 
at  which  the  potential  (or  other  function  of  position)  is  to  be  estimated 
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is  at  a  distance  from  the  nearest  portion  of  the  gravitating  mass  great 
compared  with  the  linear  dimensions  of  what  we  have  called  a  physically 
small  region.  Such  a  distance,  which,  as  we  have  seen,  may  be  extremely 
small  compared  with  the  smallest  distance  we  can  measure  directly, 
would  seem  to  mark  the  limit' of  nearness  of  P  to  the  gravitating  body 
if  the  integral  taken  for  the  hypothetical  continuous  medium  is  to  serve 
as  equivalent  to  the  true  potential.  But  further  consideration  may 
enable  us  to  push  this  limit  still  closer  to  the  body.  For  the  inac- 
curacies whose  importance  is  magnified  by  decreasing  distance  do  not, 
for  a  given  position  of  P,  occur  in  the  case  of  each  molecule  of  the 
body;  they  arise  only  in  connexion  with  the  molecules  that  are  near 
P.  Now  such  molecules,  though  perhaps  absolutely  numerous,  are 
generally  few  in  comparison  with  the  remaining  molecules  of  the  body, 
and  it  is  possible  that  their  numerical  inferiority  may  prevail  over  the 
advantage  of  their  position  in  such  a  way  as  to  render  the  total 
inaccuracy  corresponding  to  them  a  negligibly  small  fraction  of  the 
whole  potential. 

Instead  of  the  function  r~l  which  occurs  in  expressions  for  the 
potential,  let  us  consider  some  other  function  f  of  position  relative  to 
P,  which  tends,  as  r  becomes  smaller,  to  become  infinite  of  the  same 
order  as  r"^,  so  that,  for  small  values  of  r,f  is  of  the  form  kr~*  where 
k  is  a  finite  function  of  relative  angular  position.  Taking,  in  the  first 
instance,  a  single  physically  small  element  of  matter,  say  of  volume  e3, 
it  is  clear  that  the  difference  between  the  sum  ^mf  and  the  integral 
fpfdr  though  the  element  will  in  general  be  of  the  same  order  of  magni- 
tude as  either  quantity  separately  so  long  as  r  for  all  points  of  the 
element  is  of  the  same  order  of  magnitude  as  e,  but  that  the  difference 
will  diminish  to  a  quantity  smaller  in  a  ratio  comparable  with  cr"1 
when  r  becomes  great  compared  with  e.  Hence,  for  purposes  of 
estimating  order  of  magnitude,  it  is  fair  to  represent  the  difference 
between  2,mf  and  fpfdr  by  the  expression  fAer~lpfdr  where  A  is 
a  finite  number. 

To  include  all  elements  of  the  body  near  to  P,  we  suppose  the  least 
value  of  r  for  a  molecule  near  to  P  to  be  vjt  and  take  the  integral 
representing  inaccuracy  through  all  space  between  the  concentric 
spheres  r  =  ^  and  r  =  a,  where  a  is  large  compared  with  e.  If  p  is  the 
greatest  value  of  p  in  this  space,  the  order  of  magnitude  of  the  inac- 
curacy is  the  same  as  or  less  than  that  of 
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where  k'  is  a  finite  constant  replacing  k,  and  kirr^dr,  the  volume 
between  the  spheres  of  radii  r  and  r  +  dr,  takes  the  place  of  dr.  This 
is  equal  to 


of  which,  when  77  is  of  the  same  order  as  e,  the  first  term  is  small  of 
order  e  and  therefore  negligible  always,  while  the  second  term  is  small 
of  the  same  or  higher  order  provided  /x,  <  2  ;  the  second  term  would  be 
small,  but  not  of  so  high  an  order  of  smallness,  if  //.  were  between  2  and  3. 
The  case  of  ^  =  2  would  turn  on  the  order  of  magnitude  of  e  log  rj  or 
c  log  e,  which  is  very  small  though  not  as  small  as  e.  Sometimes  the 
special  form  of  the  function  k,  taken  in  connexion  with  probable  sym- 
metry in  the  average  distribution  of  molecules  round  Py  increases  still 
further  the  order  of  smallness. 

It  follows,  therefore,  that  if  ^  <  2  the  inaccuracy  is  certainly  as 
negligible  as  eft"1,  and  that  if  /*<  2  \  the  inaccuracy  is  certainly  as 
negligible  as  x/e^r1,  where  b  is  some  length  which  is  physically  not  small, 
e.g.  a  centimetre.  The  case  of  ^  =  1  corresponds  to  potential  ;  for 
attraction  components  /*  =  2. 

5.  Thus  it  appears  that  the  representation  of  potentials  and 
attractions  by  means  of  integrals  extended  through  the  hypothetical 
continuous  medium  which  replaces  the  actual  gravitating  body  is  valid 
without  sensible  error  not  only  for  points  well  outside  the  body,  but 
also  for  points  whose  distance  from  the  nearest  portion  of  the  body  is 
small  of  the  order  of  the  physically  small  length  e.  This  includes  the 
case  when  P  is  so  close  to  the  apparent  outer  surface  of  the  body  as  to 
be  sensibly  just  not  in  contact  with  it,  and  also  the  case  when  P  is  in 
a  small  but  not  imperceptibly  small  cavity  cut  in  the  body,  that  is 
a  cavity  of  such  a  size  that  the  piece  excavated  would  have  the  pro- 
perties of  matter  in  bulk  rather  than  the  properties  of  a  few  molecules. 

As  might  be  expected,  any  attempt  to  justify  the  use  of  the  same 
integral  expressions  for  the  potential  and  attractions  at  a  point  P 
which  is  at  a  distance  from  the  nearest  molecule  of  a  higher  order  of 
smallness  than  e,  results  in  failure.  For  now  a  simple  molecule 
contributes  to  the  potential  (for  example)  a  term  mr~l  which,  in  spite 
of  the  smallness  of  m,  may  become  very  great  as  r  diminishes  ;  how 
small  r  may  become  we  cannot  say,  its  least  possible  value  must 
depend  on  the  extent  to  which  the  '  impenetrability  '  of  matter  is 
true  of  isolated  molecules,  for,  since  potential  is  only  physically 
interpre  table  as  the  negative  potential  energy  per  unit  mass  of  a  particle 
(at  least  one  molecule)  at  P,  the  least  value  of  r  is  the  least  possible 
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distance  between  the  centres  of  two  molecules.  While  not  knowing 
this  least  value,  we  cannot  but  admit  the  possibility  that  a  few  terms 
of  the  type  mr~l  might  easily  become  so  important  as  to  make  the 
potential  quite  different  from  the  value  of  fpr~ldr  to  which,  as  we 
shall  see  later,  the  parts  of  the  hypothetical  continuous  distribution 
near  P  contribute  only  a  negligible  amount.  But  there  is  in  any  case, 
from  the  point  of  view  of  physics,  no  motive  for  pursuing  the  enquiry 
to  such  small  values  of  r,  for  there  are  reasons  for  supposing  that  the 
Newtonian  law  of  attraction  does  not  hold  good  at  such  distances.  In 
proving  that  the  ordinary  integral  representations  of  potential  and 
attractions  are  valid  for  distances  of  P  from  the  attracting  body  which 
are  indefinitely  small  from  the  point  of  view  of  molar  physics,  we  have 
done  all  that  would  be  required  to  justify  their  ordinary  use  in  the 
theory  of  gravitation. 

6.  One  point  requires  emphasis.     By  the  attraction  at  a  point  P 
inside  a  body  we  mean  the  attraction  (per  unit  mass)  on  a  molecule  or 
particle   at  P,  situated  in  .a   cavity  of  dimensions  which  are   only 
physically  small.     Hence  if  we  describe  a  closed  geometrical  surface  S, 
however  small,  in  a  body,  we  cannot  calculate  the  force  exerted  on  the 
part  of  the  body  inside  S  by  the  rest  of  the  body  by  use  of  the  attrac- 
tion  integrals.  ^  This   can   only  be   done  when   there   is   a  gap,  not 
smaller  than  physically  small,  between  the  attracting  and  the  attracted 
matter,  such  a  gap  as  might  be  made  by  cutting  the  body  along  the 
surface  S  and  keeping  the  fissure  open  so  that  the  opposite  sides  of  it 
are  nowhere  in  contact.     In  the  absence  of  such  a  physical  separation, 
account  must  be  taken  of  unknown  forces  between  molecules  that  are 
very  near  together.     When  the  volume  enclosed  by  $  is  not  small 
beyond  the  physical  limit   of  smallness,  such   molecules  will  all  lie 
relatively  near  the  surface  S,  and  the  forces  between  them  will  appear 
as  surface  forces  between  the  geometrically  separated  portions  of  the 
body.     In  fluids  the  extra  force  is  the  fluid  pressure,  in  solids  it  is  less 
simple. 

In  the  case  of  an  absolutely  continuous  body  there  is  nothing 
corresponding  to  the  limit  of  physical  smallness,  and  if  the  Newtonian 
law  were  supposed  to  hold  for  all  distances  however  small  there  would 
be  no  surface  forces  of  the  kind  described.  The  assumption  of  surface 
forces  in  the  ideal  case  of  continuity  is  really  a  tacit  assumption  that 
the  Newtonian  law  breaks  down  ultimately  as  r  diminishes. 

7.  We  might,  of  course,  devise  other  definite  integrals  than  those 
above   considered,   in   the   hope   of    representing   the   same   physical 
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quantities  with  possibly  greater  accuracy.  For  example  p  might  be 
obtained  by  averaging  through  a  region  smaller  than  physically  small, 
so  small  that  the  number  of  molecules  in  it  might  be  sometimes  two  or 
one  or  even  zero ;  in  this  case  fractions  of  molecules  would  become 
important,  and  the  question  would  arise  how  a  molecule  ought  to  be 
regarded  when  it  is  neither  altogether  inside  nor  altogether  outside  the 
region.  Again  we  might  reduce  the  region  of  averaging  to  the  limit 
of  mathematical  smallness  and  so  get  a  p  which  is  absolutely  zero  in 
intermolecular  space,  and  presumably  finite  and  continuous  in  the 
spaces  occupied  by  the  various  molecules.  Against  such  integrals  it  is  to 
be  urged  firstly  that  one  important  factor  of  the  function  to  be  integrated, 
namely  p,  is  inaccessible  to  experimental  measurement,  secondly  that 
even  if  p  were  known  the  integrals  would  probably  be  more  difficult  to 
evaluate  than  the  sum  Sww"1,  and  thirdly  that  the  greater  accuracy 
which  they  seem  to  possess  would  be  entirely  vitiated  by  the  probable 
failure  of  the  Newtonian  law  for  short  distances.  Moreover  a  method 
which  involves  integrating  through  the  volumes  of  individual  molecules, 
if  it  has  any  physical  significance  at  all,  implies  the  view  that  a 
molecule  is  of  the  nature  of  a  small  continuous  mass  whose  smallest 
parts  have  the  same  kind  of  properties  as  the  whole ;  this  view  is 
directly  contrary  to  modern  views  of  the  constitution  of  matter,  and 
the  mathematical  method  corresponding  to  it,  so  far  from  being  the 
best  possible  representation  of  the  facts,  must  share  all  the  defects  of 
the  method  of  summation  for  the  various  molecules. 


II.     Potentials  and  Attractions  of  accurately  continuous 

bodies. 

8.  The  potential  and  the  attraction  components  of  a  finite  body 
of  accurately  continuous  substance,  at  an  external  point  P,  are 
represented  by  volume  integrals  which,  for  ordinary  laws  of  density, 
give  rise  to  no  mathematical  difficulties.  The  subjects  of  integration 
are  finite  at  all  points  of  the  region  of  integration,  and  the  integrals  them- 
selves are  finite  and  differentiable  with  respect  to  the  coordinates  of  P 
by  the  method  known  as  'differentiation  under  the  sign  of  integration.' 
Thus  the  potential  integral,  defined  as  /pr~Wr,  justifies  its  right  to  the 
name  'potential'  by  possessing  the  property  that  its  differential 
coefficients  with  respect  to  the  coordinates  (£,  ??,  £)  of  P  are  the  attrac- 
tion integrals  of  the  type  Jp  (x  -  £)  r~3dr. 
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But  it  is  quite  another  thing  when  we  come  to  consider  the  potential 
and  attractions  at  a  point  inside  the  gravitating  body.  For  now,  for 
example,  if  we  define  the  potential  as  fpr~ldr  taken  throughout  the 
whole  body,  the  subject  of  integration  pr~l  becomes  infinite  at  the 
point  P,  a  point  in  the  volume  of  integration,  and  it  becomes  a  question 
whether  the  integral  symbol  represents  a  finite  quantity  at  all,  and,  if 
so,  whether  it  is  differentiate  and  what  are  its  differential  coefficients. 

These  troublesome  questions  might  be  avoided  by  introducing,  as 
in  the  investigation  for  a  body  of  molecular  structure,  a  cavity  within 
which  P  must  be  situated.  And,  indeed,  this  still  seems  to  be 
demanded  by  the  physical  -interpretation,  since  potential  and  attraction 
are  physically  denned  as  work  function  and  force  per  unit  mass  for  a 
hypothetical  small  mass  or  particle  at  the  point  P ;  such  particle 
cannot  be  supposed  to  occupy  space  already  occupied  by  other  matter, 
and  hence  must  be  situated  in  a  cavity  made  for  it.  But  whereas,  in 
the  case  of  molecular  structure,  there  was  suggested  from  physical 
considerations  a  limit  to  the  order  of  smallness  of  the  cavity  contem- 
plated, corresponding  in  fact  to  the  order  of  smallness  of  the  necessarily 
present  inaccuracy  in  the  mathematical  representation  adopted,  no 
such  limit  suggests  itself  in  the  case  of  continuous  bodies.  The 
retention  of  a  cavity,  of  any  definite  though  arbitrarily  chosen  order 
of  smallness,  is  not  demanded  when  there  is  no  limit  to  the  possible 
smallness  of  a  portion  of  matter,  and  would  moreover  involve  a  want 
of  precision  or  at  least  a  restriction  on  the  meaning  of  the  mathematical 
symbols  employed  which  would  considerably  discount  their  utility. 
Whereas  it  is  only  for  the  sake  of  mathematical  precision  that  the 
hypothetical  continuous  bodies  are  generally  made  the  subject  of  study 
in  preference  to  the  actual  molecular  bodies  of  which  they  are 
approximate  representations. 

9.  We  obtain  the  defmiteness  we  desire,  and,  as  will  be  seen, 
conform  at  the  same  time  to  the  conventions  and  definitions  of  Integral 
Calculus,  by  framing  new  definitions  of  the  potential  and  the  attraction 
components  at  a  point  P  (£,  77,  £),  inside  a  continuous  body.  We  first 
suppose  the  point  P  to  be  in  a  cavity,  we  then  make  the  cavity  smaller 
and  smaller,  and  define  the  limits  (if  such  exist)  to  which  the  potential 
and  the  attraction  components  at  P  tend  with  the  vanishing  of  the 
cavity  as  the  potential  and  the  attraction  components  respectively  at 
P  when  no  cavity  exists.  It  must  be  recognised  that  this  passage  to 
the  limit  entirely  destroys  the  physical  meaning  which  the  quantities 
considered  possess  at  any  stage  short  of  the  limit,  but  on  the  other 
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hand  it  gives  us  extremely  convenient  standard  approximations  to 
these  quantities  in  cases  of  physical  interest ;  the  very  definition  of 
the  term  limit  implies  that  the  approximation  can  be  made  as  close  as 
we  please  by  taking  the  cavity  sufficiently  small. 

It  is  also  to  be  noticed  that  a  relation  such  as  X  =  ^  (where  X 

°$ 

is  a  force  component  and  F  the  potential),  which  holds  inside  a  cavity 
of  finite  size  however  small,  might  not  persist  after  passage  to  the 
limit.  That  is  to  say,  though  of  necessity 

T-        8F 

Lim  X  =  Lim  -^  , 

9c 
it  is  not  equally  inevitable  that 

Lim  X  —  T-T  Lim  V. 
ft 

In  fact  if,  as  is  customary,  we  drop  the  phrase  'limit'  from  our 
notation,  though  keeping  the  idea  in  mind,  we  have  to  face  the  fact 

that  the  formula  X  -  -TT-  ,  valid  for  free  space,  requires  examination 

before  we  can  be  sure  that  it  is  true  at  a  point  in  the  substance 
of  the  body.  And  if  it  be  objected  that  the  formula  is  known  to 
be  true  in  all  cases  of  physical  interest,  and  that  no  such  interest 
attaches  to  its  validity  or  otherwise  in  the  case  which  has  avowedly 
no  physical  significance,  an  answer  is  that  if  we  decide  to  use  a 
certain  kind  of  mathematical  functions  as  approximate  representations 
of  physical  quantities,  we  must  become  acquainted  with  the  meanings 
and  properties  of  these  functions  before  we  can  make  intelligent  use 
of  them. 

Hence  it  is  natural  for  the  student  of  the  theory  of  attractions 
to  turn  his  attention  to  that  part  of  pure  mathematics  which  has  to 
do  with  the  definition  and  properties  of  volume  and  surface  integrals. 


III.     Volume  integrals. 

10.  Let  /  be  a  function  of  position,  and  let  a  finite  volume  T  be 
divided  into  a  great  number  of  elements  AT,  of  small  linear  dimensions  ; 
let  yi  be  a  quantity  associated  with  an  element  of  volume  AT,  chosen 
according  to  some  law,  so  that  it  is  either  the  value  of  f  at  some  point 
of  the  element,  or  at  any  rate  not  greater  than  the  greatest  or  less  than 
the  least  value  of  f  for  points  in  the  element.  If  f  is  finite  at  all 
points  in  the  volume  T,  the  sum  2/iAr  extended  to  all  elements  of  T 
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is  finite,  and  will  remain  so  no  matter  how  small  and  correspondingly 
numerous  are  the  elements  AT.  If  this  sum  tends  to  a  limit  as  the 
number  of  elements  tends  to  infinity,  and  the  linear  dimensions  of  each 
tend  to  zero,  arid  if  this  limit  is  independent  of  the  law  specifying  / 
and  of  the  manner  of  subdivision  into  elements,  the  limit  is  called  the 
volume  integral  of  /  through  the  volume  T,  and  is  denoted  by  ffdr. 
This  definition  is  only  valid  on  the  supposition  that  /  is  finite  at  all 
points  in  T. 

Whether  the  limit  here  spoken  of  does  or  does  not  exist  depends  on 
the  nature  of  the  function/;  we  shall  assume  that  it  does  exist  for  all 
the  forms  of/  which  we  meet  with  in  potential  theory. 

11.  Next  consider  the  case  in  which  /  is  a  function  which  becomes 
infinite  at  a  point  P  within  the  volume  T]  clearly  we  need  a  new 
definition,  and  that  which  has  been  generally  adopted  is  as  follows. 
Surround  the  point  P  by  a  small  closed  surface  t,  and  take  the  volume 
integral  through  the  whole  of  the  volume  T  except  the  part  included 
by  t ;  we  thus  exclude  P  from  the  range  of  integration,  and  so  get  a 
finite  integral.  Now  let  the  surface  t  become  smaller  and  smaller, 
whilst  always  surrounding  P  ;  if  the  volume  integral  tends  to  a  finite 
limit  as  the  space  enclosed  by  t  tends  to  vanishing,  and  if  the  limit  is 
independent  of  the  shape  of  t,  then  this  limit  is  defined  to  be  the 
integral  of  /  throughout  the  whole  volume  T.  The  definition  may  be 
expressed  symbolically  thus  : 


[T  rT 

I   fdr  =  Lim  I   fdr, 

J  t^oJt 


where  the  symbol  -*•  is  used  to  denote  such  phrases  as  'tending 
towards'  or  'tends  towards/  so  that  t  -*-  0  reads  'as  t  tends  towards 
zero.'  Here  and  elsewhere  the  subscript  to  the  integral  specifies  the 
inner  boundary  of  the  region  of  integration. 

If  we  call  the  space  inside  the  vanishing  surface  t  a  'cavity'  in 
the  volume  of  integration,  we  see  at  once  the  parallelism  between  the 
definition  of  this  kind  of  volume  integral  and  that  of  the  so-called 
potential  and  attractions  at  a  point  in  the  substance  of  a  continuous 
body. 

The  volume  integral  (if  it  exists)  through  a  region  within  which  / 
becomes  infinite  at  some  point  is  seen,  by  the  above  definition,  to  be  a 
mathematical  conception  of  a  different  character  from  the  integral  for 
a  region  in  which  /  is  everywhere  finite.  In  a  sense  one  might  say 
that  the  latter  is  a  true  volume  integral  while  the  former  is  the  limit 
of  a  true  volume  integral.  The  latter  bears  to  the  former  the  kind  of 
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relation  that  a  single  limit  bears  to  a  double  limit,  or  that  a  finite 
series  bears  to  the  so-called  sum  of  an  infinite  series. 

12.  Analogously  with  the  terminology  of  series,  we  speak  of  the 
volume  integral  as  convergent  if  it  tends  to  a  finite  limit  with  the 
vanishing  of  the  cavity,  divergent  if  it  tends  to  become  infinitely  great, 
and  semi-convergent  if,  as  sometimes  happens,  there  is  a  finite  limit 
whose  value  is  not  independent  of  the  shape  of  the  cavity.     Divergent 
integrals  are,  for  ordinary  purposes,  as  meaningless  as  divergent  infinite 
series,  and  so  we  must  satisfy  ourselves  that  the  integrals  in  use  in 
gravitation  problems  are  convergent  either  absolutely  or  in  the  con- 
ditional manner  corresponding  to  semi-convergence. 

To  decide  whether,  for  a  given  form  of  /,  the  integral  is  convergent 
or  not,  we  have  the  following  rule,  depending  on  the  order  of  the 
infinity  of  /  at  P  in  terms  of  r,  the  distance  from  P  to  the  point  at 
which  f  is  estimated.  If  /  becomes  infinite  at  P  of  an  order  lower 
than  r~3  the  volume  integral  is  convergent,  if  of  an  order  higher  than 
r~3  the  integral  is  divergent ;  if  of  the  order  r~9  exactly  the  integral 
may  be  divergent,  semi-convergent,  or  convergent,  according  to  the 
way  in  which  f  in  the  neighbourhood  of  P  depends  on  the  angular 
position  of  r.  This  rule  is  not  stated  with  sufficient  accuracy  to  rank 
as  a  theorem,  and  one  can  easily  think  of  exceptions  to  it;  for  example 
the  case  of/=  r~4  cos  0  (in  the  notation  of  spherical  polar  coordinates), 
which  is  obviously  convergent  for  any  cavity  symmetrical  about  the 
plane  0  =  \tr  though  otherwise  likely  to  be  divergent,  shews  that  some- 
thing like  semi-convergence  may  be  associated  with  infinities  of  order 
greater  than  3 ;  but  the  rule  is  a  convenient  approximation  to  the  facts. 

13.  With  a  view  to  justifying  the  rule  here  given,  it  will  be 
convenient  to  re-state,  with  slight  modification  of  form,  the  definition 
of  convergence.      The  integral  of  /  through  the  volume   T,  which 
includes  a  point  of  infinity  P,  is  convergent  if,  corresponding  to  any 
arbitrarily  chosen  small  quantity  cr,  there  can  always  be  found  a  closed 
surface  0  surrounding  P  such  that  all  closed  surfaces  t  surrounding  P 
and  lying  wholly  inside  6  have  the  property  that 


fdr 


<  (T. 


That  this  is  essentially  the  same  as  the  definition  of  §  11  appears  at 
once  when  we  think  of  the  ordinary  definition  of  a  limit ;  for*  if  the 
integral  through  T  has  a  limit  A,  we  can  choose  0  so  that 


<  * 
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and  therefore 


<cr; 


it  is,  of  course,  to  be  understood  that  P  must  not  lie  on  the  surface  0. 
Thus  the  property  constituting  the  definition  of  convergence  of  the 
present  Article  is  a  consequence  of  the  property  laid  down  as  a 
definition  in  §  11. 

Conversely,  possession  by  an  integral  of  the  property  specified  in 
the  present  Article  involves  as  a  necessary  consequence  the  existence  of 
a  limit  A,  though  giving  no  indication  of  its  actual  value.  For  by 
taking  0  sufficiently  small  we  can  keep  the  fluctuation  of  the  value  of 
the  integral  for  different  cavities  within  6  as  small  as  we  please,  that  is 
small  without  limit  ;  and  infinitely  restricted  fluctuation  is  the  same  as 
infinite  approximation  to  some  definite  (and  therefore  finite)  value. 

Part  of  the  rule  of  §  12  maybe  formulated  in  the  following  theorem. 
If  within  a  sphere  of  finite  radius  (a),  having  P  as  centre,  f  is  everywhere 
less  in  absolute  value  than  Mr~^t  where  M  is  a  (finite)  constant  and 
/x<3,  the  integral  is  convergent.  To  prove  this  let  us  take  for  the 
surface  0  the  sphere  r  =  r),  where  'n<a)  and  let  us  denote  by  «  the 
distance  from  P  to  the  nearest  point  of  the  surface  t  of  the  cavity  ;  the 
cavity  is  of  course  entirely  inside  0,  but  is  otherwise  unrestricted  as  to 
shape.  Since  the  modulus  of  a  sum  is  not  greater  than  the  sum  of  the 
moduli,  and  since  an  integral  is  the  limit  of  a  sum, 


M 


where  the  subscript  e  means  that  the  inner  boundary  of  the  integrals 
is  the  sphere  r  —  e ;  the  second  inequality  holds  because  \f\  is  positive 
and  e  is  completely  inside  0. 

In  dealing  with  a  function  of  r  only,  we  may  combine  all  elements 
dr  that  lie  between  spheres  of  radii  r  and  r  +  dr  in  the  single  expression 
so  that 


\f. 


fdr 


3-f* 
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it  being  noted  that  e  <  77  and  that  3  -  /x  is  positive,  so  that  the  last 

expression  obtained  is  positive.     Now  if  a-  be  any  arbitrarily  chosen 

i 

small  quantity,  we  have  only  to  take  r)  less  than  {(3  -  ft)  a/lvM}*-*  in 
order  to  get  a  surface  0  such  that 

I  C° 

fdr    «r, 
\J  t 

whatever  shape  t  may  have  provided  only  it  lies  inside  0.  Thus  the 
convergence  of  the  integral  of /is  established. 

It  will  be  noticed  that  the  convergence  of  the  integral  of  /  in 
accordance  with  this  theorem  involves  also,  as  the  proof  indicates,  the 
convergence  of  the  integral  of  |/|. 

It  need  hardly  be  pointed  out  that  the  position  and  shape  of  the 
outer  boundary  T  of  the  region  of  integration  do  not,  in  general,  aifect 
the  question  of  convergence  ;  whatever  the  outer  boundary  may  be, 
provided  it  does  not  include  other  points  of  infinity,  it  is  only  the 
part  of  the  volume  just  round  P  that  is  in  danger  of  making  the 
integral  very  great,  and  so  only  that  part  need  be  studied  with  a  view 
to  detecting  divergence. 

It  is  clear  that  the  theorem  holds  equally  well  for  cases  in  which 
the  point  P  where  the  infinity  occurs  is  not  inside  but  just  on  the 
boundary  of  the  region  of  integration. 

14.  The  corresponding  theorem  for  divergence  is  as  follows.  If 
within  a  sphere  of  finite  radius  (a),  having  P  as  centre,  f  is  everywhere 
algebraically  greater  than  mr~*  where  m  is  a  constant  greater  than  zero, 
and  /x  ^  3,  the  integral  is  divergent.  To  prove  this,  we  take  as  outer 
boundary  the  sphere  r  -  a,  and  as  inner  boundary  a  surface  t,  and  we 
denote  by  e  the  distance  from  P  to  the  furthest  point  of  the  surface  t, 
so  that  the  sphere  r  =  e  completely  surrounds  the  cavity.  Then 

ra  fa  ra 

I  fdr  >  m   I    r~^  dr>  m  \   r~^dr  • 
Jt  J  t  Je 

as  before,  we  collect  all  the  elements  dr  between  r  and  r  +  dr  into  the 
expression  kirr^dr,  and  so  get 

I    fdr  >  ^Trm    I    r*~^  dr, 

Jt  7c 

;-(/*- 3)  _a-(M-3)}    or    47rm  {log  »  -  log  e}, 


according  as  /x  is  greater  than  or  equal  to  3.  In  either  case  the 
expression  obtained  tends  to  infinity  for  e  -^-  0 ;  and  so  the  integral, 
being  greater  than  a  quantity  which  tends  to  become  endlessly  great, 
is  divergent. 
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15.  The  case  of  semi-convergence  need  only  be  illustrated  by  an 
example.     Suppose  that  f  is  r~3  cos#,  and  consider  the  values  of  the 
integral  in  the  regions  having  a  common  outer  boundary  r  =  a,  and 
having  for  inner  boundaries  in  the  first  instance  the  sphere  r  =  e,  and 
in  the  second  instance  the  sphere  r2  -  re  cos  0  =  2€2,  these  being  two 
small  spheres  of  which  the  latter  touches  and  completely  surrounds  the 
former. 

The  difference  between  the  integrals  over  these  two  regions  is  the 
integral  through  the  space  between  the  two  small  spheres,  which  is 
certainly  not  zero  so  long  as  c  is  different  from  zero,  since,  if  we  consider 
the  volumes  of  the  region  cut  off  by  a  cone  of  small  solid  angle  having 
the  origin  as  vertex,  the  positive  contribution  to  the  integral  from  the 
frustum  where  cos  0  is  positive  is  greater  in  absolute  value  than 
the  negative  contribution  from  the  frustum  where  cos  0  is  negative. 
Further,  the  magnitude  of  the  integral  is  independent  of  e,  since  if  we 
multiply  e  by  k  we  can  get  the  new  region  of  integration  by  multiplying 
all  radii  vectores  from  P  by  k,  and  thus  each  element  of  volume  dr  is 
multiplied  by  P ;  the  subject  of  integration  r~3  cos  6  is  correspondingly 
multiplied  by  k~3,  and  so  the  integral  is  unaltered.  Thus  the  integral 
over  the  space  between  the  small  spheres,  being  finite  when  €  is  not 
zero,  has  the  same  finite  value  as  e  tends  to  vanishing  ;  in  other  words 
there  is  a  finite  difference  between  the  values  of  the  original  integral 
corresponding  to  the  different  cavities.  It  is  clear,  from  the  symmetry 
of/  about  the  plane  0  =  JTT,  that  the  integral  is  zero  for  the  cavity  whose 
centre  is  P,  and  therefore  not  zero  for  the  cavity  whose  centre  is  not 
at  P  ;  but  in  neither  case  is  it  infinite.  Thus  the  semi-convergence  of 
the  integral  is  demonstrated. 

This  example  suggests  the  remark  that  two  cavities  are  to  be 
regarded  as  of  different  shapes  even  if  they  are  similar,  if  they  are  not 
similarly  situated  with  respect  to  P.  The  cavities  considered  in  the 
example  are  both  spheres,  but  since  one  has  P  at  its  centre  while  in 
the  other  P  trisects  a  diameter,  the  cavities  are  regarded  as  of  different 
shapes  for  purposes  of  the  present  discussion. 

IV.     Theorems  connecting  volume  and  surface  integrals. 

16.  There  is  a  well-known  theorem  connecting  volume  integrals 
with   surface   integrals   taken   over   the   boundary   of    the   region   of 
volume-integration.     If  the  region  be  finite,   if  /,  m,  n  denote  the 
direction  cosines  of  the  normal  drawn  outwards  from  the  region  at 
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a  point  of  the  boundary  B,  and  if  £,  r/,  £  be  functions  having  finite 
space  differential  coefficients  at  all  points  in  the  region, 

............  (ft 


where  c?$  represents  an  element  of  area  of  the  boundary,  the  surface 
integral  is  taken  over  the  complete  boundary,  and  the  volume  integral 
through  the  complete  volume.  A  proof  is  given  in  Williamson's 
Integral  Calculus,  Chapter  XL 

It  is  worth  while  to  enquire  whether  this  theorem  can  be  extended 
to  the  case  in  which  there  is  a  point  P  in  the  volume  where  the 
subject  of  volume-integration  becomes  infinite.  The  course  which 
suggests  itself  is  to  surround  the  point  P  by  a  small  closed  surface  o-, 
and  to  apply  the  original  theorem  to  the  region  bounded  internally 
by  a-  and  externally  by  the  surface  B.  The  complete  boundary 
consists  of  both  B  and  <r,  and  so  we  get 


where  the  suffixes  to  the  surface  integrals  indicate  the  surfaces  over 
which  they  are  taken.  Now  if  the  subject  of  integration  of  the 
volume  integral  is  such  as  to  make  it  convergent  with  respect  to  the 
infinity  at  P,  the  left-hand  side  of  the  equality  tends  to  a  definite 
limit  as  the  dimensions  of  o-  decrease  towards  zero.  Consequently  we 
get  as  the  limiting  form  of  the  equality, 


=  {  (l£  +  mri  +  n£)dS+Um  (  (It  +  mrj  +  nQdS  ...(2). 

J  B  (T-^OJ  <r 

When  the  volume  integral  is  convergent,  the  left  side  of  (2)  is  a 
perfectly  definite  finite  quantity,  and  hence  the  limit  indicated  on  the 
right-hand  side  must  exist  and  be  independent  of  the  shape  of  o-  ;  it 
will  exist,  but  have  a  value  dependent  on  the  shape  of  o-,  if  the  volume 
integral  is  semi-convergent.  In  the  former  case  it  is  frequently  con- 
venient to  determine  the  value  of  the  limit  by  taking  some  specially 
simple  shape  for  <r,  such  as  a  sphere  with  centre  at  P.  Generally,  if 
the  subject  of  volume  integration  is,  in  the  neighbourhood  of  P,  of 
order  P~**(/n<8),  where  r  denotes  distance  from  P,  the  subject  of  the 
surface  integral  under  the  limit  sign  is  of  order  r~^+1  where  r  equals 
the  radius  e  of  the  sphere  o-  ;  also  dS  is  e2c?w,  where  du>  is  an  element 
of  solid  angle,  and  so  the  surface  integral  is  of  order  €3~"  and  tends  to 
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the  limit  zero  for  e  ^-0.  If  ^  =  3,  the  case  of  possible  semi-convergence, 
the  surface  integral  is  of  order  e°,  so  far  as  its  dependence  on  e  is 
concerned,  and  therefore  may  have  for  limit  a  value  different  from 
zero. 

17.  A  generalisation,  and  at  the  same  time  a  particular  case,  of 
the  fundamental  '  surface  and  volume  integral  theorem '  is  got  by 
putting  $£,  <f>r),  cf>£,  instead  of  £,  ->?,  £,  where  <j>  is  another  function  of 
position  which  has  finite  space  differential  coefficients  at  all  points  in 
the  region.  The  volume  integral  then  becomes 


so  that  the  theorem  takes  the  form 


Now  <£  is  continuous  and  therefore  finite  throughout  the  whole 
region  of  integration,  tut  let  us  suppose   that  some   or  all  of  the 

functions  £,  17,  £,  —  ,    r-  become  infinite  at  a  point  P  in  the 

QX     oy     oz 

region.  If  the  infinities  are  such  that  both  the  volume  integrals  are 
convergent,  it  is  clear  that  by  introducing  the  cavity  a  and  making  it 
tend  to  zero  dimensions  we  get  the  relation 


I    </>.(/£  +  my  +  n£)  dS  +  Lim  I    <£  .  (/£  +  mrj  +  n£)  dS 

JB  a-^OJo- 


[ 

=1 


and  that  when  the  convergence  is  due  to  £,  y,  £  being  infinite  of  lower 
order  than  r~2  the  limit  of  the  surface  integral  is  zero.  When  the 
volume  integrals  are  semi-convergent,  or  when  £,  >?,  £  are  infinite  of 
the  same  order  as  r~2,  the  limit  of  the  surface  integral  may  be 
different  from  zero,  possibly  depending  on  the  shape  of  the  cavity  ; 
it  would  usually  be  convenient  to  take  it  in  the  ultimately  equivalent 
form 

(5), 


where  <£p  signifies  the  value  of  <£  at  the  point  P, 

2—2 
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18.     Green's  Theorem  is  got  from  the  *  surface  and  volume  integral 
theorem  '  by  putting 

dV  dV  dV 


where  U,  V  are  functions  of  position  ;  if  we  use  the  notation  r-  for 

differentiation  along  the  outward  normal,  and  A  for  Laplace's  operator 

a2      a2      9 


whence 


(6) 


by  symmetry.  These  two  equalities  constitute  Green's  Theorem.  The 
statement  of  the  theorem  requires  modification  when  the  region  of 
integration  is  multiply  connected  and  either  U  or  F  is  many-  valued  ; 
this  modification  is  discussed  in  Maxwell's  Electricity  and  in  Lamb's 
Hydrodynamics,  and  need  not  be  entered  upon  here. 

If  one  of  the  functions,  say  U,  together  with  all  of  its  differential 
coefficients  which  occur  in  the  formula  (6),  is  finite  throughout  the 
region,  and  if  F  becomes  infinite  at  a  point  P  in  the  region,  we 
isolate  P  in  a  cavity  a-  and  make  the  dimensions  of  o-  tend  to  zero.  If 
the  volume  integrals  converge,  we  get  a  pair  of  equalities  precisely 
similar  to  the  above,  save  that  to  the  first  member  we  must  add 

Lim  I   U  -^-  dS,  and  to  the  third  member  Lim  I  V-^-dS. 

v^Qjo-  OV  tr-^OJo-         & 

Generally  the  convergence  of  all  the  volume  integrals  would  involve 
that  V  should  become  infinite  at  P  of  an  order  lower  than  r~l,  since 
when  F  is  a  function  of  position  relative  to  P,  as  is  frequently  the 
case,  a  space  differentiation  adds  one  to  the  order  of  the  infinity  so 
that  A  F  is  of  an  order  higher  by  r~2  than  F.  In  this  case  both  the 
integrals  over  the  surface  a-  are  of  the  order  of  a  positive  power  of  the 
small  length  r,  and  their  limits  are  zero. 

The  case  of  V-  r~l  is  one  of  special  interest  in  physical  applications; 
it  is  also  interesting  mathematically  because  it  is  just  the  case  in 
which  semi-convergence  is  to  be  looked  for.  There  is  no  semi- 
convergence  however,  for  jU^(r~^dr  is  absolutely  zero,  since 
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A  (r~1)=0  at  all  points  between  a-  and  the  outer  boundary,  and  the  other 
volume  integrals  are  convergent  because  the  subjects  of  integration  are 
infinite  of  lower  order  than  r~s  ;  and  if  all  the  terms  but  one  of  an 
equality  are  definite  or  convergent,  that  one  cannot  be  semi-convergent. 
Clearly,  since  dS  is  comparable  with  r*d<a,  where  do>  is  an  element  of 

solid  angle,  /    r~l  -r-  dS  has  a  zero  limit  ;  but  /    U-r-  (r~l)  c?$(when  a- 
J  &        w  J  o     w 

is  a  sphere  of  radius  e,  which,  in  the  absence  of  semi-convergence, 
may  be  assumed  without  loss  of  generality)  has  the  same  limit  as 


or  - 


that  is   Upjdu  or 

Thus  Green's  Theorem  in  this  case  gives  the  equalities 


(7). 


Almost  identical  reasoning  applies  to  the  case  in  which  V  satisfies 
AF=0  at  all  points  of  the  region  other  than  P,  and  becomes  infinite 
at  P  in  such  a  way  that  Lim  (r  F)  =  Mt  where  M  is  a  definite  constant. 

r-*0 

The  theorem  becomes 


Another  case  of  Green's  Theorem  much  used  in  Physics  is  that  in 
which  U=V.     The  two  equalities  reduce  to  the  single  one 


no  particular  interest  attaches  to  an  examination  of  possible  modifica- 
tions in  this  formula  when  V  has  an  infinity  at  a  point  in  the  region 
of  integration. 

V.    The  differentiation  of  volume  integrals. 

19.  We  shall  next  discuss  the  possibility  of  differentiating  a 
volume  integral  with  respect  to  a  parameter  which  occurs  in  the  subject 
of  integration  but  does  not  affect  the  boundary  of  the  region  of 
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integration.  The  only  parameters  that  need  be  considered  here  are 
the  coordinates  of  a  point  P  at  which  the  subject  of  integration  /  has 
an  infinity.  We  shall  call  these  coordinates  £,  77,  £,  keeping  at,  y,  z  to 
denote  the  coordinates  of  the  element  dr  of  integration.  The  question 
to  be  settled  is  whether  the  integral  has  a  differential  coefficient  with 
respect  to  £,  and  if  so  whether  that  differential  coefficient  is  equal  to 
the  integral  of  dffi£.  Integration  means  passage  to  a  limit,  and  so  also 
does  differentiation ;  we  have  to  find  out  whether  alteration  of  the 
order  of  the  two  passages  to  limit  alters  the  value  of  the  final  result. 
When  P  is  in  the  region  of  integration  there  is  in  each  case  the 
additional  passage  to  limit  corresponding  to  the  closing  of  the  cavity, 
and  the  question  to  be  settled  is  whether  by  first  integrating,  second 
closing  the  cavities,  third  making  A  £  ->  0,  (where  A  £  is  an  increment 
of  £),  we  get  the  same  result  as  by  first  making  A  £  -*-  0,  second  in- 
tegrating, third  closing  the  cavity. 

20.  First  we  shall  consider  the  case  in  which  P  is  outside  the 
region  of  integration,  and  shew  that  if /has  at  all  points  of  the  region 
and  for  all  contemplated  values  of  £  a  differential  coefficient  with 
respect  to  £,  which  differential  coefficient  is  a  uniformly  continuous* 
function  of  £  throughout  the  region,  the  differential  coefficient  of  the 
integral  is  the  same  as  the  integral  of  the  differential  coefficient. 

The  incremental  ratio  of  the  integral  is 


/ 


{/(£+**)- 

which,  by  the  theorem  of  mean  value, 


where  *  =/'  (£  +  0  A  £)  -/'  (£),  and  1  >  0  >  0. 

Now  the  uniform  continuity  of  f'(£)  implies  that  for  an  arbitrarily 
chosen  small  quantity  or  we  can  always  find  a  quantity  w  such  that  for 
all  values  of  A  £  less  than  w, 


and  consequently  also  |  e  ,  is  less  than  o-,  for  all  points  in  the  region  T  of 
integration.  Thus  by  choosing  A  £  less  than  o>  we  can  ensure  that  |  fcdr  \ 
shall  be  less  than  o-  T,  which,  in  virtue  of  the  finiteness  of  77and  the  arbi- 
trariness of  o-,  is  arbitrarily  small.  In  fact  the  difference  between  the 

*  It  can  be  proved  that  if  a  function  is  continuous  at  all  points  in  a  region  it  is 
uniformly  continuous  throughout  the  region. 
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integral  of  /'(£)  and  the  incremental  ratio  of  the  integral  of  /(£)  can 
be  made  arbitrarily  small.    Hence  the  limit  of  the  incremental  ratio,  i.e. 

~7  //(£)  dr,  is  equal  to  //'(£)  dr,  as  we  set  out  to  prove. 

21.  Next  we  consider  the  case  in  which  P  is  within  the  region  of 
integration.  The  rough  rule  for  this  case  is  that  if  the  original  integral 
is  convergent,  and  if  the  integral  obtained  by  differentiating  under  the 
sign  of  integration  is  convergent,  the  latter  is  the  differential  coefficient 
of  the  former. 

It  will  serve  our  purpose  to  prove  this  proposition  for  a  particular 
case,  namely  that  in  which  the  subject  of  integration  /  is  the  product 
of  two  factors,  each  subject  to  special  restrictions.  One  of  these,  which 
we  denote  by  p  (#,  y,  z)  or  briefly  by  p,  is  supposed  to  be  a  function  of 
absolute  position,  not  involving  £,  17,  £  at  all  ;  it  is  assumed  to  be  finite, 
not  to  vanish  at  P,  and  to  have  space  differential  coefficients  which  are 
uniformly  continuous  throughout  the  region  of  integration.  The  other 
factor  is  supposed  to  be  a  function  of  position  relative  to  P,  and  to 
become  infinite  at  P  ;  we  may  denote  it  by  </>  (x  —  £,  y  -  77,  z  —  £),  or 
sometimes  for  brevity  by  <J>  (£,  x]  or  <£  (£).  It  has  obviously  the 
property  that  <£  (£  +  A  £,  x)  -  $  (£,  x  —  A  £),  so  that 


rT 

The  integral  to  be  differentiated  is  I    p.<f>.dr,  or,  written  in  full, 

rT 
Lim  I    p.<f>.dT,  where  e  is  a  cavity  surrounding  the  point  P.      The 

e-9-OJe 

incremental  ratio  is 


i  r       ( 

--Lim 

A£|_e'-H».0./e' 


where  e'  is  a  cavity  surrounding  the  point  P'  (£  +  A  £  y,  £),  which  may 
be  taken  to  be  in  all  respects  similar  to  e.  To  get  the  differential 
coefficient  of  the  integral  it  is  necessary,  in  the  incremental  ratio,  first 
to  make  e  and  e'  tend  to  vanishing,  and  afterwards  to  make  A  £  -»-  0. 

Before  we  pass  to  either  limit,  however,  we  are  at  liberty  to  simplify 
the  form  of  the  incremental  ratio  in  any  manner  that  seems  desirable. 

/T 
P  (x,  y,  z)  .  <f>  (£  +  A  £)  .  dr  imagine  the  boundaries  of 
e 

the  region  of  integration,  and  every  volume  element,  shifted  a  distance 
A  £  in  the  negative  direction  of  the  axis  of  x.  An  element  of  volume 
originally  at  a  point  K'  is  merely  shifted  to  a  point  K  whose  position 
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relative  to  the  point  (£,  17,  £)  is  the  same  as  that  of  K'  relative  to 
(£  +  A  I,  17,  0  ;  so  that  if  .AT  is  (#,  y,  2),  the  value  of  <f>  belonging  to  the 
shifted  volume  element,  originally  presenting  itself  as  <f>  (a?+  A  £,  £+  A  |), 
is  equally  well  represented  by  the  form  <£  (#,  £)  ;  but  the  value  of  p 
for  the  element  now  brought  to  K  is  that  appropriate  to  K\  i.e. 
p  (a?  +  A  £,  ^,  z).  The  inner  boundary  e'  of  the  integral  is  brought  by 
the  shifting  into  coincidence  with  e,  but  the  outer  boundary  is  changed 
to  a  surface  T'  which  is  simply  T  displaced  without  change  of  form. 
In  fact 

rT  rT' 

I  rj»(ai).+(6+*t9a)dr=j    P(®+  A  £).<£(£,#)  dr, 

so  that  the  incremental  ratio  of  the  integral  with  as  yet  unclosed 
cavity  is  equal  to 


or  to 


where  the  latter  integral  is  extended  to  the  region  between  T  and  T', 
being  taken  positively  where  the  boundary  of  T'  lies  outside  T, 
negatively  where  the  reverse  is  the  case.  By  the  theorem  of  mean 
value  the  above  expression  equals 


where  w  =  px'  (x  +  0  A  £,  y,  2;)  -  p^  (#,  y,  0),  and  1  >  6  >  0. 

Now  as  P  is  not  on  the  boundary  of  T,  A  £  can  always  be  taken 
small  enough  to  prevent  P  from  being  in  the  region  between  Tand  T'  ; 
hence  the  integral  for  this  region  has  no  dependence  on  the  cavity  c. 
The  assumed  uniform  continuity  (which  includes  fmiteness)  of  px't  and 
the  assumed  convergence  of  the  original  integral,  are  sufficient  guaran- 
tees of  the  convergence  of  the  integrals  of  px'  </>  and  o><£.  Hence  we 
may  proceed  to  the  limits  for  e-^-O  and  e'-^O,  and  so  get  the  relation  : 

FT 

Incremental  ratio  of    /    p^>  dr 

A  £).<£(£).  dr..  .(9), 


and  the  cavities  e,  e'  are  now  closed  up  and  finished  with. 
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As  A  £  becomes  smaller  it  is  clear  that  the  volume  between  T  and 
T  approximates  to  a  very  thin  shell  over  the  surface  of  T  whose  normal 
thickness  (outwards  from  T)  is  -  A  £ .  /,  where  /  is  the  x  cosine  of  the 
outward  normal.  Thus  the  corresponding  volume  integral  approxi- 
mates to  and  has  the  same  limit  as  the  surface  integral 


-Aff 

h 

And  in  virtue  of  the  uniform  continuity  of  px'}  corresponding  to  any 
arbitrary  small  quantity  o-,  we  can  always  find  a  quantity  K  such  that 
for  all  values  of  A  £  less  than  K,  and  for  all  points  of  the  region  of  inte- 
gration, to  |  <  o-,  and  so 

I   rT  rT 

\     o><f>dr    «r  I     \<j>\dr. 

This  last  expression  is  o-  multiplied  by  a  finite  quantity,  for,  since  p  is 
finite  and  does  not  vanish  at  P,  the  convergence  of  /  \<j>\dr  is  implied 

rT 

in  the  assumed  convergence   of   I    p<f>dr,  at  least  if  the  latter  be 

rT 

convergence  of  the  kind  discussed  in  §  13.     Hence    I   <o<£  dr  can,  by 

suitable  choice  of  A  £,  be  made  smaller  than  any  arbitrary  small 
quantity,  and  so  tends  to  the  limit  zero  for  A  £  -*•  0.  Proceeding  now 
to  the  limit  A  £  -^-  0  in  relation  (9),  we  get 

(10). 


The  theorem  of  §  17,  formula  (4),  enables  us  to  transform  the  right- 
hand  side  of  (10),  giving 


9  [T 
—  j 


provided  the  last  integral  is  convergent. 

If  the  function  p  is  of  such  a  simple  character  near  P  that  the 
nature  of  the  integrals  depends  entirely  on  the  form  of  <f>,  it  is  clear 
that  the  case  of  possible  semi-  convergence  is  covered  by  the  above 
reasoning,  certainly  as  far  as  formula  (10);  but  in  formula  (11)  the 
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use  of  formula  (4)  may  introduce  an  extra  term  on  the  right-hand 
side,  namely  the  limit  of  the  surface  integral  of  lp$  over  a  new  cavity 
round  P.  One  can  imagine  that  peculiarities  in  the  form  of  />  might 
invalidate  some  of  the  steps  of  the  argument,  but  such  peculiarities  are 
not  to  be  expected  in  physical  applications. 

VI.     Applications  to  Potential  Theory. 

22.  The  potential  at  a  point  P  (£,  17,  £)  of  a  finite  mass  of 
continuous  matter,  whose  density  at  a  point  (#,  y,  z)  is  p,  a  function 
of  xy  y,  z  but  not  of  £,  r/,  £,  is  the  volume  integral  V  =  fpr~ldr,  where 


r  =  J^  (x  ~  £)2  ;  the  attraction  component  parallel  to  the  axis  of  x  is  X 
where  X  =  j(x  -  fjfr~3  dr  ;  both  integrals  are  taken  through  the  whole 
space  occupied  by  the  body. 

When  (£,  •»;,  £)  is  outside  the  body,  and  p  is  finite  and  subject  to 
such  restrictions  as  are  required  for  the  validity  of  the  theorems  proved 
in  the  preceding  Articles,  there  is  no  infinity  of  the  subjects  of  inte- 
gration in  the  region,  and  so 


a    3 
Also  = 


so  that 

32  v     22  y     ww       r     /  32        ^2        32  \ 

V-V\  fJ-r  — 


23.  When  the  point  P  is  inside  the  body,  the  potential  and  the 
attraction  components  have  no  longer  the  simple  physical  interpreta- 
tion suggested  by  their  names,  but  are  defined  as  the  limits  of  these 
physical  quantities  in  a  vanishing  cavity.  And  this,  as  we  saw  in 
§  11,  implies  their  equivalence  to  the  integrals  represented  by  the 
same  symbols  as  in  §  22,  but  referring  to  a  region  including  the  point 
of  infinity  P  of  the  subjects  of  integration,  and  therefore  only 
intelligible  when  the  integrals  are  convergent. 

The  subject  of  the  potential  integral,  being  infinite  at  P  of  the 
order  of  r~l,  the  integral  is  convergent  ;  and  the  attraction  integrals, 
having  subjects  of  integration  that  are  infinite  of  the  order  of  r~2,  are 
also  convergent.  Hence  we  have,  by  §  21, 
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But  if  we  differentiate  X  with  respect  to  £  under  the  sign  of 
integration,  we  get  an  integral  whose  subject  of  integration  is  of  the 
order  of  r~3,  so  that  there  is  a  possibility  of  semi-convergence  or 
divergence.  Instead,  therefore,  of  merely  quoting  a  simple  differen- 
tiation rule  in  this  case,  we  must  proceed  with  care. 

It  is  clear  that  for  the  integral  X  =  fp  (x  —  £)  r~sdr  the  argument 
of  §  21  holds  as  far  as  the  formula  (10),  which  in  this  case  becomes 


the  volume  integrals  involved  being  convergent  and  all  cavities  being 

o  -y 

closed  up.     This  formula  shews  that  -j-?  has  a  definite  value. 

Now  surround  P  by  a  small  surface  a-  and  use  formula  (4)  of  §  17, 
putting  p  for  the  quantity  there  called  </>,  and  (x  -  £)  r~3  for  the 
quantity  there  called  £.  Thus  we  get 


pLim  I  l(x- 

(T-*.0  J   <T 


whence 

o  -p-  ,.  /.y        3 

±  -  Pp  Urn  l-l(*-.Q  r~3dS  -  Lim  /    pf  {(x  -  £)  r~s]  dr  .  .  .(13). 

d£  <r-*»OJ<r  <r-*~oJ<r        W 

The  sum  of  the  limits  here  indicated  is  perfectly  definite  and  inde- 
pendent of  the  shape  of  o-,  but  either  limit  taken  separately  has  a  value 
which  depends  on  the  shape  of  a-  ;  this  can  be  seen  readily  by  studying 
the  surface  integral  first  when  o-  is  a  sphere  and  second  when  o-  is  a 
very  flat  cylinder  with  plane  ends  parallel  to  the  plane  of  x. 

Since  formulae  corresponding  to  (13)  hold  for  I7'  and  ^/,  we  get 
by  addition 

(Li         9  Jc         S^y 


Now  here  the  subject  of  volume  integration  is  identically  zero  at  all 
points  outside  the  cavity,  and  so  the  integral  is  zero  whatever  the 
shape  of  the  cavity,  and  its  limit  is  zero.  Hence  the  value  of  the 
surface  integral  is  independent  of  the  shape  of  the  cavity,  and  may  be 
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calculated  on  the  assumption  that  or  is  the  sphere  r  =  e  ;  in  this  case 
/=-(#-£)  e-1,  so  that  SJ  (#-£)  =  -  e,  and  dS=  €2du>,  where  du  is  an 
element  of  solid  angle  at  P  ;  thus  the  integral  becomes  -  /<&>,  which 
equals  -  4?r.  Thus  our  equality  becomes 

ax    8F    az 


a2F    82F    82F 
W  +  W  +  W  =    ^pp' 
which  is  Poisson's  equation. 

24.  The  theorems  proved  above  for  the  differential  coefficients  of 
F  are  perfectly  intelligible  for  a  body  of  the  hypothetical  continuous 
structure  which  we  have  postulated..  But  when  applied  to  a  body  of 

molecular  structure  such  a  symbol  as  -p-  requires  qualification.    It  has 

been  seen  that  for  a  continuous  body  A  £  was  only  made  to  tend  to  zero 
after  we  had  first  closed  the  cavities  e  and  e'  corresponding  to  F  and 
F  +  A  F  ;  in  fact  A  £  though  small,  was  always  large  compared  with 
the  dimensions  of  the  cavities,  but  this  fact  did  not  interfere  with  our 
making  A  £  as  near  to  zero  as  we  pleased. 

But  for  a  body  of  molecular  structure  the  cavities  must  always  be 
large  enough  to  be  capable  of  containing  a  great  number  of  molecules, 
and  so  we  can  never  close  them  entirely  ;  hence  A  £,  so  far  from  ever 
vanishing,  must  be  actually  large  compared  with  the  smallest  length 
which  can  be  regarded  as  only  physically  small  ;  nevertheless  A  £  may 
be,  to  our  senses,  extremely  small.  Hence  instead  of  the  true 

differential  coefficient  -53-  we  have  the  incremental  ratio  —  ~-r  .  where 
of  A£ 

A  £  though  very  small  is  still  definitely  prevented  from  attaining  the 
higher  orders  of  smallness  which  lie  on  the  way  to  the  limit  zero. 

3F 
Thus  it  is  clear  that  the  symbol  -^  ,  just  as  V  itself,  is  inexact  and 

stands  for  something  not  precisely  defined  ;  but  the  inaccuracy  or 
deviation  from  a  precise  value  is  no  greater  than  the  inaccuracy  which 
regards  matter  as  continuous,  and  is  in  fact  an  inaccuracy  so  small  as 
to  be  inappreciable  to  our  senses.  Accordingly  the  relations 

X—  -££    and    3  ~-r  —  —  47TD 

9c  8£ 

have  a  sufficiently  precise  meaning  when  applied  to  bodies  of  molecular 
structure. 
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VII.    Applications  to  Theory  of  Magnetism. 

25.  If  a  body  is  magnetised  so  that  the  components  of  the 
intensity  of  magnetisation  at  a  point  (#,  y,  z)  are  A,  B,  C,  the 
magnetic  potential  at  an  external  point  P,  (£,  rj,  £),  is  given  by 


and  the  x  component  of  magnetic  force  is  a  where 


the  body  being  regarded  as  of  continuous  structure  ;  so  long  as  P  is 
outside  the  body  it  is  clear  that 


Outside  the  body  the  induction  (a,  b,  c)  is  the  same  as  the  force 
(a,  (3,  y),  and  therefore  remembering  that  A,  B,  C  are  functions  of 
set  y,  z  but  not  of  £,  rj,  £,  while  r  depends  only  on  x-£,y-v],  2  -  £,  we 
see  that 

a  =  -  ((A  £-  +  B  ~  +  C^-\  ^  (r-1)  dr 


dy2          dz2          cxty         dx 

=  [(0*  -  A  I-}  I  (r->)  dr-ffAJ—S  pi  j-  (r->)  dr 
J\dx  dz/dz^  )\      fy          W  Sy  ^      ' 

—A%-)  I  (r->)  dr  +  ((A  ~-B  1}  I-  (r-)  d 

dx  dz/  3C  J\      ty  dxj  8rj  ^ 


-  Sr, 

where 


jP,  G,  H  are  the  components  of  the  vector-potential  at  P  ;  the  relation 
between  induction  and  vector  potential  is  frequently  written 

,G,ff)  .....................  (19). 
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26.  When  P  is  inside  the  magnetised  body  the  integral  of 
formula  (14)  is  convergent,  and  so  the  formula  may  stand  as  the 
definition  of  the  potential  at  P.  The  properties  of  V,  thus  defined, 
are  most  easily  deduced  from  another  expression,  obtained  by  making 
a  cavity  round  P  and  applying  the  theorem  of  §  17,  formula  (4),  taking 
0  to  be  r~l  and  writing  A,  B,  C  for  £,  77,  £.  The  surface  integral  over 
the  cavity  has  a  zero  limit,  and  so  we  get 


where  the  surface  integral  refers  only  to  the  boundary  T  of  the  body, 
and  the  cavity  is  now  closed  and  finished  with.  This  form  of  the 
magnetic  potential  exhibits  it  as  equivalent  to  the  gravitation 
potential  of  a  volume  distribution  of  density 

_M_a^?_!<7 

p        da      ty      dz  ' 

combined  with  a  distribution  of  surface  density  I  A  +  mB  +  nC  spread 
over  the  boundary  of  the  body.  (Of  course  formula  (20)  is  equally 
true  when  P  is  outside  the  body.) 

If  we  suppose  that  P  is  right  inside  the  body,  i.e.  not  on  the 
boundary,  there  is  no  infinity  in  the  subject  of  surface  integration  ; 
the  volume-integral  part  of  V  has  the  properties  of  the  gravitation 
potential  studied  in  Section  VI.  Thus  V  has  definite  space  differential 
coefficients  obtained  by  differentiating  under  the  sign  of  integration 
in  formula  (20)  (not  in  formula  (14)*);  accordingly,  since 


a  (T/dA     dB     dC\  a 

nC)  —  (r-1)  dS-  /    [  ^-  +  ^~  +  -5-  )  5-  (r~ 

'  dx^  /     \fa       dy        dzjdx^ 


the  volume  integral  having  a  perfectly  definite  value.  Now  we  cut 
a  cavity  o-  round  P,  and  apply  the  theorem  of  §  17,  formula  (4),  and 
we  get 

8F  f   .  8 

— TT  —  —  Lim  I  (IA  +  mB  +  nC)  r-  (r 

8%  <r-*-OJ<T  V& 

+  Lim  I    ( A  ^-  +  B  —  - 

o-^-O  Jo-  \       3®  ty 

These  two  limits  combined  give  a  value  which  is  independent  of  the 
shape  of  o-,  but  the  value  of  each  limit  taken  separately  depends  on  the 

*  The  theorem  of  §  21  does  not  apply  to  the  integral  of  formula  (14),  since  the 
infinity  at  P  is  of  the  order  r~2. 
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shape  of  o-;  the  volume  integral  we  see,  by  comparison  with  (15),  to  be 
the  x  component  of  the  magnetic  force  in  the  cavity  due  to  all  the  matter 

outside  the  cavity.  So  in  general  —  -^  is  not  the  limit  of  the  com- 
ponent of  force  in  the  cavity,  but  differs  from  it  by  an  amount  repre- 
sented by  the  limit  of  the  surface  integral.  If,  however,  we  can  choose 
a  shape  for  the  cavity  which  shall  make  the  limit  of  the  surface  integral 
zero,  the  limit  of  the  force  component  in  the  cavity  will  be  accurately 

represented   by  -  ^7-  ;    and  this   is   effected  by   making  the   cavity 

a  cylinder  whose  generators  are  parallel  to  the  direction  of  the  vector 
(A,  B,  (7)  at  the  point  P,  with  flat  ends  perpendicular  to  the 
generators,  all  the  linear  dimensions  of  the  cylinder  tending  to  zero 
in  such  fashion  that  the  linear  dimensions  of  the  ends  tend  to  become 
vanishingly  small  compared  with  the  length ;  this  may,  for  brevity,  be 
called  a  'long'  cylinder.  The  direction  chosen  for  the  generators 
ensures  that  the  integral  of  I  A  +  mB  +  nC  for  the  curved  portion  of 
the  surface  tends  to  zero,  and  the  relative  smallness  of  the  flat  ends 
makes  the  integral  over  these  tend  also  to  zero.  The  definition  of  the 
magnetic  force  (a,  /?,  y)  at  a  point  P  in  the  body  is  '  the  limit  of  the 
force  in  a  cavity  in  the  form  of  a  long  cylinder  with  generators  parallel 
to  the  resultant  intensity  of  magnetisation';  and  this  definition,  in 
connexion  with  the  present  argument,  justifies  the  statement  that 

dV 

nr 

The  definition  of  the  induction  (a,  b,  c)  at  a  point  P  in  the  body  is 
'  the  limit  of  the  force  in  a  cavity  in  the  form  of  a  very  flat  circular 
cylinder  with  generators  parallel  to  the  resultant  intensity  of  magnetisa- 
tion,' where  by  a  very  flat  cylinder  is  meant  one  whose  linear  dimensions 
tend  to  zero  in  such  a  way  that  the  length  tends  to  become  vanishingly 
small  in  comparison  with  the  linear  dimensions  of  the  plane  ends.  For 
such  a  cavity  IA  +  mB  +  nC  tends  to  zero  on  the  curved  part  of  the 
surface,  to  -  /  over  one  of  the  plane  ends,  and  to  +  /  over  the  other, 
/  being  the  resultant  intensity  of  magnetisation ;  and  each  of  these  ends 
ultimately  subtends  a  solid  angle  2?r  at  P.  Thus  the  three  surface 
integrals  of  which  that  in  (21)  is  a  type  have  for  limits  the  components 
of  force  at  a  point  between  two  infinite  circular*  parallel  planes,  the 

*  The  word  '  circular '  is  introduced  in  order  to  exclude  cases  in  which  the 
resultant  force  at  a  point  between  the  parallel  planes  is  not  normal  to  them.  The 
circles  are  supposed  to  have  a  common  axis,  passing  through  P. 
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one  covered  with  a  uniform  surface  density  /,  the  other  with  a  uniform 
surface  density  —  /,  of  matter  that  attracts  according  to  the  Newtonian 
law  ;  this  force  is  known  to  be  47r/  perpendicular  to  the  planes,  and  so 
its  components  are  47rJ.,  farB,  473-  (7.  So,  for  the  flat  cavity,  (21)  yields 
the  equality 

a  =  -4:TrA  +  a  ...........................  (22). 

27.  The  integrals  of  formulae  (18)  representing  vector  potential 
are  convergent  for  a  point  inside  the  body,  and  may  therefore  stand 
as  the  definition  of  the  vector  potential  at  such  a  point.  If  in  the 
formula  (4)  of  §  17  we  put  0  for  £,  -  C  for  77,  B  for  £,  and  r~l  for 
<£,  we  get 

F=jT(nB-mVrW-f(^-%)r-><h  ......  (23), 

the  cavity  of  formula  (4)  being  closed  and  finished  with,  and  the 
surface  integral  over  the  cavity  having  a  zero  limit.  This  result 
exhibits  F  as  the  gravitation  potential  of  a  finite  volume  distribution 
combined  with  a  surface  distribution  ;  it  shews,  therefore,  that  F  has 
definite  differential  coefficients  with  respect  to  the  coordinates  of  P. 
From  the  two  formulae  analogous  to  (23), 

dff      dG         f  \  I  m        t     \   d         IA  r>7^    8 

-=  --  -&  =  /  \(Cl 
8>7        8£      jTl 


[T\fiA      8£\   8       /dC     BA\  81 

(  -  --  ^—  )  =  --  (  5  --  •*—  )  5-     r     dr.  .  .(24). 
J    |_\83/       9#  )  dy      \dx       dzjdz] 


Cut  a  cavity  o-  round  P  and  apply  the  theorem  of  formula  (4)  §  17, 
and  the  result  is  readily  seen  to  be 

-  Bl         r~*  dS 


=  _  Lim  f  \(Cl  -An)?--  (Am-  Bl)  f\ 
,^-0J«L  3*  fyj 


wherein  the  limits  on  the  right-hand  side  together  give  a  value 
independent  of  the  shape  of  o-,  though  the  value  of  each  separately 
depends  on  the  shape  of  cr.  The  volume  integral  is  the  same  as 


and  accordingly  represents  the  x  component  of  force  due  to  all  the 
magnetisation  outside  the  cavity ;  the  surface  integral,  together  with 
the  corresponding  surface  integrals  in  the  two  other  formulae  analogous 
to  (25),  will  tend  to  zero  if  A/l  =  Bjm  =  C/n  over  practically  the  whole 


26-28]  SURFACE   INTEGRALS  33 

surface  of  the  cavity,  and  this  is  ultimately  the  case  when  the  cavity 
is  the  flat  cylinder  used  in  defining  the  induction.  For  this  shape  of 
cavity  (25)  is  equivalent  to 

dH    dG 


which,  with  the  two  other  equalities  of  the  same  type,  constitutes  the 
vector  relation 

(a,  b,  c)  =  curl  (Ft  G,  H)t 

true  now  for  points  inside  as  well  as  for  points  outside  the  magnetised 
body. 

It  should  be  noticed  that  the  definition  of  vector  potential  used 
in  the  present  discussion  is  not  that  which  is  regarded  as  fundamental 
in  the  physical  theory,  though  equivalent  to  it.  The  usual  definition 
is  contained  in  the  relation  (19)  coupled  with  the  relation 


It  is  easy  to  verify,  on  the  lines  of  the  present  Article,  that  the 
vector  defined  by  the  relation  (18),  whether  for  internal  or  for  external 
points,  satisfies  this  further  condition. 

VIII.     Surface  Integrals. 

28.  When  gravitating  matter  is  distributed  in  a  very  thin  layer, 
or  when  the  surface  of  a  body  is  charged  with  electricity,  the  corre- 
sponding potential  and  attraction  at  a  point  P  are  represented  by  surface 
integrals,  a  surface  density  a-  taking  the  place  of  a  volume  densit}7". 
The  integrals  are  of  the  type  fafdS  where  a-  is  usually  free  from  such 
mathematical  peculiarities  as  might  raise  doubts  concerning  the 
existence  of  the  integrals,  and  /  is  a  function  having  an  infinity  at  the 
point  P  (£  ??,  0- 

So   long   as   P  is  not  actually  in  the  region  of  integration  the 

integrals  do  not  present  any  difficulties,  and  the  formulae  X  =  -^t 


=  0,  are  clearly  valid. 

When  the  point  P  is  in  the  surface  distribution  we  must  cut 
a  cavity  round  it  of  dimensions  that  tend  to  zero,  and  the  question 
of  convergence  necessarily  arises.  We  consider  first  the  case  in 
which  the  integration  takes  place  over  a  portion  of  a  plane  surface. 

The  chief  test  of  convergence  is  now  as  follows.  If  within  a  circle 
of  finite  radius  (a),  having  the  point  P  as  centre,  the  subject  <£  of 
L.  3 
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integration  is  always  less  in  absolute  value  than  Mr~^  where  /A  <  2  and 
M  is  a  definite  constant,  the  integral  f<J>dS  is  convergent.  To  prove 
this  we  shall  shew  that,  corresponding  to  any  arbitrarily  chosen  small 
quantity  oy  there  can  always  be  found  a  closed  curve  0  surrounding  P 
such  that  all  closed  curves  t  surrounding  P  and  lying  wholly  inside 
0  have  the  property  that 


Take  for  the  curve  0  the  circle  r  =  17,  where  rj  <  a,  and  denote  by  e  the 
distance  from  P  to  the  nearest  point  of  the  boundary  t  of  the  cavity; 
the  cavity  is  of  course  entirely  inside  0,  but  is  otherwise  unrestricted 
as  to  shape.  Since  the  modulus  of  a  sum  is  not  greater  than  the  sum 
of  the  moduli, 


dS, 


,  . 

2  —/A  2  —/A 

1 

Hence  by  choosing  rj  less  than  {(2  —  p)  <r/27rM}  '2~»  we  get  a  carve 
0  satisfying  the  specified  condition  ;  the  integral  is  accordingly 
convergent. 

When  the  order  of  the  infinity  of  <jf>  is  the  same  as  that  of  r~2, 
semi-convergence  may  appear. 

29.  Passing  to  the  case  in  which  the  region  of  integration  is  a 
portion  of  a  curved  surface,  we  shall  assume  P  to  be  a  point  at  which 
there  is  a  definite  tangent  plane  and  such  that  at  all  points  of  the 
region  within  a  finite  distance  of  P  the  principal  curvatures  are  both 
finite.  We  need  consider  only  the  integral  taken  through  a  finite 
region  not  extending  far  from  P,  and  in  virtue  of  the  finiteness  of  the 
curvatures  at  and  near  P  we  can  always  choose  this  region  so  that, 
if  Q  is  any  point  of  it  and  0  the  inclination  of  the  tangent  plane  at  Q 
to  the  tangent  plane  at  P,  for  all  positions  of  Q  in  the  region  0  <  a, 
where  a  is  a  definite  acute  angle.  Let  the  projection  of  Q  on  the  tangent 
plane  at  P  be  Q0,  let  r,  r0  represent  PQ,  PQ0  respectively,  dS  an 
element  of  area  round  Q,  dS0  the  projection  of  dS  on  the  tangent  plane 
at  P,  B  the  boundary  of  the  area  of  integration,  and  JB0  its  projection 
on  the  tangent  plane  at  P. 
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Since  dS  =  d8Q  sec  0, 

B  fBn 

<£  sec  0  dS0) 


the  second  integral  being  taken  in  the  tangent  plane  at  P. 
If  within  the  region  of  integration 

\  +  \<M9r* 
where  M  is  a  constant  and  2  >  /*  >  0,  then 

|  0  sec  6  |  <  Mr-*  sec  0 


IS 


or,  since  r0<r,  and  sec  0  < sec  a, 

|  <£  sec  0 1  <  M  sec  a  r0~*, 

where  M  sec  a  is  finite  since  a  is  acute. 

/•.BO 
Hence   I     <£  sec  0  dS0    is    convergent,    and    therefore    so    also 

rS 

I    <j>d8',  thus  the  test  of  convergence  is  the  same  whether  the  surface 

of  integration  be  plane  or  curved  provided  the  curvatures  be  finite. 
The  existence  of  a  definite '  tangent  plane  at  P  is  not  a  necessary 
feature  in  the  proof,  the  essential  thing  is  that  there  shall  be  a  finite 
region  round  P  for  which  0  <  a  <  JTT,  6  being  inclination  to  some  fixed 
plane  through  P\  for  example  the  surface  might  be  a  cone  and  P  its 
vertex.  (Compare  Poincare',  Potentiel  Newtonien,  §  33.) 

30.  Applying  the  test  of  the  preceding  Articles  we  see  that  at 
a  point  in  a  surface  distribution  of  gravitating  matter  or  electricity  the 
potential  is  represented  by  a  convergent  integral,  but  the  attraction 
components  in  the  tangent  plane  are  represented  by  integrals  whose 
order  renders  semi-convergence  possible.  It  is  not  difficult  to  shew, 
by  a  particular  example,  that  semi-convergence  does  occur ;  for  the 
attraction  of  a  uniform  plane  elliptic  disc  (of  eccentricity  e  and  surface 
density  a-)  at  a  focus  is  27r<r(l  -  \/l  -e*)/e  if  the  cavity  is  circular,  but 
is  zero  if  the  cavity  is  an  ellipse  similar  and  similarly  situated  to  the 
edge  of  the  disc,  with  the  focus  for  centre  of  similitude  ;  the  verification 
of  these  statements,  by  using  polar  coordinates  and  integrating,  is 
quite  easy. 

The  component  of  attraction  at  P  normal  to  the  surface  is 
represented  by  a  convergent  integral,  but  this  quantity  is  the  attraction 
in  a  cavity,  though  a  vanishing  one,  and  must  be  distinguished  from 
the  normal  component  of  attraction  at  a  point  very  close  to  the 
unbroken  surface  but  not  in  it ;  it  is,  in  electrical  applications,  the 

3—2 
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'mechanical  force  per  unit  charge/  the  quantity  denoted  by  Rz  in 
Prof.  J.  J.  Thomson's  Elements  of  Electricity  and  Magnetism,  §  37, 
whereas  the  normal  attraction  at  a  point  just  not  in  the  surface  is  the 
quantity  there  denoted  by  R. 

31.  The  distinction  drawn  above,  between  the  attraction  at 
a  point  in  the  surface  and  that  at  a  point  just  not  in  the  surface,  brings 
us  to  a  question  of  a  kind  which,  for  lack  of  a  fourth  dimension, 
does  not  arise  geometrically  in  the  case  of  volume  integrals,  the 
question,  namely,  whether  an  integral  /<£  dS  tends  to  a  definite  limit 
if  the  point  P,  where  <£  has  an  infinity,  is  not  originally  in  the  surface, 
but  approaches  a  point  of  the  surface  as  a  limiting  position. 

Let  0  be  the  point  of  the  surface  to  which  P  gets  continually 
nearer ;  it  will  be  convenient  to  take  0  as  origin  of  coordinates  and  the 
tangent  plane  at  0  as  plane  of  z ;  we  shall  suppose  that  there  is  a 
limiting  position  of  the  line  PO,  as  P  moves  up  to  coincidence  with  0, 
which  makes  with  the  plane  of  z  a  definite  angle  different  from  zero 
and  so  has  definite  direction  cosines  lot  mot  n0,  of  which  the  last  is 
numerically  greater  than  zero.  The  length  PO  will  be  denoted  by 
K,  and  the  coordinates  of  P  by  (£,  17,  £)  or  (-1*,  —niK,  —m<\  while 
#,  y,  z  represent  the  coordinates  of  a  variable  point  Q  on  the  surface ; 
the  subject  of  integration,  <£  (x,  y,  z,  £,  rj,  £),  may  for  brevity  be 
represented  by  <£,  while  <f>  (x,  y,  z,  0,  0,  0),  the  value  of  <f>  when  P 
is  coincident  with  0,  will  be  represented  by  <£0.  If  integration  be 
extended  to  a  finite  part  of  the  surface  round  0,  bounded  by  a  closed 
curve  B,  the  quantities  to  whose  different  meanings  and  possibly 
different  values  it  is  desired  to  draw  attention  are  respectively 


rB 

I    fa 

J 


and  Lim 


The  first  thing  to  notice  is  that,  while  the  integral  of  <£  requires  no 
cavity  so  long  as  K  is  different  from  zero,  which  is  the  case  at  all  stages 
in  the  passage  to  limit  denoted  by  K  -^  0,  the  integral  of  <£0  is  only 
intelligible  in  terms  of  a  cavity  €  round  the  point  0,  though  this  cavity 
of  course  tends  to  vanishing.  If,  therefore,  we  set  out  to  find  the 
algebraic  difference  between  the  two  quantities  which  form  the  subject 
of  discussion  (which  may  conveniently  be  denoted  by  D)  we  have 


rB  rB 

DEE  Lim       <f>d8-       <f>0dS 
K-*~oJ  J 

rB  rB 

=  Lim  I     <£</$- Lim  I    <f>0d/S. 

K-*-OJ  e^-0./e 
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Since  the  term  involving  the  limit  for  K-^O  in  no  way  depends  upon 
c,  and  the  term  involving  the  limit  for  e  -^-  0  in  no  way  depends  upon 
K,  it  is  a  matter  of  indifference  in  what  order  we  suppose  the  passages 
to  limit  to  be  made  ;  accordingly  we  are  at  liberty,  if  we  please,  to 
make  first  the  passage  to  the  limit  for  K  -^  0,  so  that  at  any  stage  short 
of  the  limits  we  shall  think  of  K  as  extremely  small  compared  with  the 
linear  dimensions  of  the  cavity  e.  The  difference,  then,  before  passage 
to  either  limit,  may  be  put  in  the  form 

fids*  f*<l>dS-f  fad8. 

Now  if  we  proceed  first  to  the  limit  for  K  -->•  0,  the  points  P  and  0  at 
all  stages  of  this  passage  are  quite  outside  the  region  of  integration  of 
the  last  two  integrals,  and  the  functions  <£  and  <£0  are  kept  definitely 
removed  from  their  infinite  values  ;  hence  in  the  absence  of  peculiarities 
of  </>  other  than  that  infinity  at  P  which  is  the  special  subject  of  our 

rB 

investigation,  we  get  the  same  limit  for    I    <j>dS  whether   we  first 

integrate  and  then  make  K  -*•  0  or  first  make  K  -*•  0  and  then  integrate. 
In  fact 

Lim  I   <j>dS  =  (   Um<l>d8=  (    fadS  ; 

K-^-O  Je  Je     K-*-0  Je 

whence 

D=lAm[lAm  f  <£<£#+  Lim  /  <f>dS-  I 

e-*0  LK-S-O./  K^QJe  Je 


-Lim  Lim 


im  f'<j>d8  .......................................  (27), 

-0  J 


the  notation  implying  that  e  is  kept  constant  while  K  -*-0,  thus  yielding 
a  limit  which  is  a  function  of  e,  and  that  afterwards  the  limit  of  this 
function  of  e  is  taken  for  e^-0. 

Let  us  suppose  <£  to  be  of  the  form  (z  -  £)A  r~*  where  r  denotes  PQ 
and  A  and  /*  are  positive,  and  let  us  proceed  to  make  a  closer  examina- 
tion of  D  for  this  particular  case.  We  shall  assume  that  the  principal 
curvatures  of  the  surface  are  finite  at  all  points  in  a  finite  region  round 
0,  and  that  the  cavity  e  is  determined  by  the  intersection  of  the  surface 
with  the  narrow  cylinder  a?  +  y2  =  *2',  and  we  picture  to  ourselves 
a  small  piece  of  the  surface,  which  we  may  call  the  '  cap,'  bounded  by 
this  curve  whose  projection  on  the  plane  of  z  is  a  circle  of  radius  e  and 
centre  0,  and  a  point  P  at  a  distance  K  from  0  which  is  extremely  small 
compared  with  e.  To  begin  with,  we  observe  that  there  is  a  finite 
constant  a  such  that  for  all  points  Q  in  the  cap  \z\<  as2,  where  s  stands 
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for  v#2  +  #2,  the  distance  of  Q  from  the  axis  of  z  ;  for,  since  the  surface 
has  finite  curvature  at  all  points  of  the  cap,  |  z\/s*  tends  to  a  finite  limit 
for  any  given  azimuth  as  Q  approaches  0,  and  so  is  finite  at  all  points 
of  the  cap  ;  and  the  various  values,  being  finite,  have  a  finite  superior 
limit  a  which  is  of  the  same  order  of  magnitude  as  the  greatest 
curvature  of  a  normal  section  through  0  ;  thus  the  inequality  is 
proved. 

Consider  now  the  curve  of  intersection  of  the  surface  with  the 
cylinder  #2  +  y2  ^  /c2  ;  this  divides  the  cap  into  two  regions,  an  inner 
region  whose  linear  dimensions  are  of  the  order  of  K  and  therefore 
small  in  comparison  with  those  of  the  cap,  and  an  outer  region 
comprising  most  of  the  cap,  in  which  there  is  no  point  whose  distance 
from  0  is  of  a  higher  order  of  smallness  than  K.  The  integral  whose 
limit  is  D  can  be  regarded  as  the  sum  of  two  integrals,  one  over  the 
inner  region,  one  over  the  outer  region,  and  these  will  be  considered 
separately. 

Taking  first  the  outer  region,  and  remembering  the  assumption 
w0^0,  we  see  that,  while  there  may  be  points  in  this  region  for  which 
r  <  s,  there  are  no  points  in  it  for  which  s/r  becomes  infinite,  so  that 
there  is  a  finite  superior  limit  (3  for  the  values  of  s/r  in  the  region  ; 
the  finiteness  of  the  curvature  is  a  guarantee  that  there  is  a  finite 
superior  limit  c  (differing  from  unity  by  a  quantity  of  the  order  of  e'2) 
to  the  secant  of  the  inclination  to  the  axis  of  z  of  the  normal  to  the 
surface  at  Q,  i.e.  the  ratio  of  dS  to  its  projection  dS0  on  the  plane  of  z  ; 
and  |  £  \/s  is  finite  at  all  points  of  the  region  and  has  therefore  a  finite 
superior  limit  y',  so  that  |  £  |  <  y's  ;  as  |  z  \  <  as2,  it  follows  that 

|  Z  -  £  |  S~l  <  y  +  as  <  y 
where  y  is  a  finite  quantity. 
Thus  in  the  outer  region 


dS<cdS0t 
and  so 


I  f 


+dS 


The  limit  of  this  for  K  -^  0  and  e  -»•  0  is  zero  provided  /x  -  X  <  2. 
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Taking  next  the  inner  region,  we  know  that  in  it  z  <  as2 
while  £  is  of  the  same  order  of  smallness  as  K,  so  that  |  z  -  £  |  is  of  the 
same  order  as  K,  and  there  must  be  an  inequality  |  z  -  £  |  <  g<  where  g  is 
a  definite  constant.  Clearly  there  is  a  superior  limit  to  the  secant  of  the 
angle  between  the  normal  at  Q  and  the  line  PQ,  so  that  there  is  an 
inequality  dS  <  hr^du,  where  d<a  represents  an  element  of  solid  angle  at 
P  ;  and  r  always  bears  to  K  a  finite  ratio,  so  that  there  is  a  double 
inequality  pK>r>  <?K,  where  p  and  q  are  constants.  Hence  in  the 
inner  region 

((z  -  £)A  r-*  dS  <  gxq  ~  ^p2h^  -  »+*  (^ 

and  the  limit  of  this,  for  /c  -*•  0,  is  zero  provided  /x,  —  X  <  2. 

Thus  D  =  0  for  <£  =  (2  -£)*/•->*  subject  to  /A-A.<2,  and  it  is  an 
immediate  inference  that  D  =  Q  for  <f>  =  o-(z-£)^r~'J-  subject  to  the 
same  condition,  if  o-  is  a  function  of  #?,  y,  z  which  is  finite  throughout 
the  region  of  integration.  For  this  type  of  integral  the  condition  for  the 
vanishing  of  D  is  not  the  same  as  the  condition  for  the  convergence  of 
the  integral  of  <£0,  for,  in  the  neighbourhood  of  0,  z  becomes  small  of 
the  order  of  r02,  so  that  the  condition  of  convergence  of  <f>0  is 

//.  -  2A  <  2. 

Powers  of  x  —  £,  y  —  r]  might  appear  in  <£  ;  they  would  count  as  the 
corresponding  powers  of  r  in  applying  the  test  just  proved,  though  of 
course  they  would  not  be  equivalent  to  powers  of  r  if  one  were 
examining  a  case  where  something  analogous  to  semi-convergence 
seemed  probable. 

For  the  potential  integral  X  =--  0,  ^  =  1,  and  therefore  D  =  Q.  Thus 
the  potential  at  0  is  the  same  as  the  limit  of  the  potential  at  P  as 
P  approaches  0,  so  that  V  is  not  discontinuous  at  points  on  the 
surface. 

32.  For  the  attraction  components  /x-  —  A  =  2,  and  the  test  of  the 
previous  Article  is  not  applicable  so  that  a  special  investigation  is 
required.  Let  us  consider  first  a  tangential  component,  say  that 
whose  subject  of  integration  is  a-  (x  -  £)  r~'s. 

The  corresponding  integral  of  <£0  is  semi-convergent,  and  it  may 
appear  useless  to  investigate  the  difference  D  between  the  unknown 
limit  of  the  integral  of  <£  and  the  quantity  of  uncertain  value  which  is 
the  integral  of  c£0.  But  the  integral  of  <f>  is  quite  independent  of  the 
shape  of  the  cavity,  in  fact  it  does  not  require  any  cavity,  so  we  are  as 
free  as  in  the  case  of  absolute  convergence  to  give  to  the  cavity  any 
shape  we  please,  so  long  as  we  attach  to  the  integral  of  </>0  the  value 
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associated  with  that  particular  shape  ;  thus  the  semi-convergence  does 
not  introduce  any  uncertainty  into  the  meaning  of  D. 

Let  us  take  the  same  cavity  as  in  the  preceding  Article,  using  the 
same  notation  and  applying  whatever  parts  of  the  reasoning  remain 
valid  for  the  changed  form  of  <£.  And  let  us  consider  what  error  is 
introduced  into  the  subject  of  integration  if  we  replace  o-  by  o-0,  its 
value  at  0,  dS  by  its  projection  dS0  on  the  plane  of  z,  and  r  by  /  the 
distance  from  P  to  the  projection  of  Q  on  the  plane  of  z.  The  error 
due  to  the  change  in  a-  corresponds  to  the  omission  of  a  factor  which 
differs  from  unity  by  a  quantity  of  the  order  of  smallness  of  s,  if  we 
assume  the  function  a-  to  have  no  troublesome  peculiarities*  at  0  ;  the 
error  due  to  the  change  in  dS  corresponds  to  the  omission  of  a  factor 
differing  from  unity  by  a  quantity  of  the  order  of  s2 ;  and,  since 
|  r  -  r'  <  z\<  as2,  the  error  due  to  the  change  in  r  corresponds  to 
the  omission  of  a  factor  differing  from  unity  by  a  quantity  of  the 
order  sV"1.  So  the  most  important  terms  representing  error  in  the 
integral  over  the  cap  are  of  the  order  of  the  integrals  over  the  cap  of 
(x  -  £)  r~ss  and  (x  -  £)  r~V ;  for  these  error  integrals  X  =  0,  //-  =  1,  s  in 
the  numerator  counting  as  equivalent  to  r  since  there  is  a  finite 
superior  limit  to  sr~l,  and  so,  by  the  previous  Article,  the  error  has  a 
zero  limit  for  K  -H*.  0  and  e  ^  0. 

Hence,  for  the  x  component  of  attraction, 

D  =  Lim  Lim  I  <TO  (x  -  £)  r'~3dS0 , 

e-^O  K^0  J 

where  it  is  clear  that  the  integral  is  now  taken  over  a  circular  area  of 
radius  e  in  the  plane  of  z,  and  represents  the  component  of  attraction 
at  P  of  a  circular  disc  of  uniform  surface  density  o-0.  Now,  before 
passage  to  the  limit,  P  is  not  in  such  a  position  of  symmetry  that  the 
x  attraction  component  must  vanish,  but  if  we  describe  the  reflexion 
with  respect  to  the  plane  x  =  £  of  the  lesser  of  the  two  arcs  into  which 
this  plane  divides  the  circumference  of  the  disc,  we  obtain  a  division 
of  the  disc  into  two  areas  the  greater  of  which,  on  account  of  the 
symmetry  of  the  position  of  P  with  respect  to  it,  contributes  nothing 
to  the  attraction  component.  The  component  is  therefore  that  due  to 
the  crescent-shaped  smaller  area,  whose  mass  is  very  nearly  4£eo-0  and 
whose  nearest  point  is  at  a  distance  from  P  comparable  with  € ;  thus  the 

*  If  further  precision  be  desired,  we  may  assume  j  a-  -  <r0  \  <  Msm,  where  M  is 
finite  and  m  positive.  The  error  corresponding  to  this  is  less  than  the  integral  of 
Msm  (x  -  £)  r~3,  for  which  X  =  0,  fjt,  =  2-m,  and  the  limit  of  the  error  is  zero.  In  the 
text  m  is  taken  to  be  unity. 
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component  of  attraction  is  of  the  same  order  of  magnitude  as  <r0£e~l, 
which  has  a  zero  limit  if  we  make  *  -*•  0  (i.e.  £  -*•  0)  before  e  -»-  0.  Thus 
D  is  zero,  so  that  the  x  component  of  attraction  of  the  whole  surface 
at  P  tends  to  a  limit,  as  P  approaches  0,  equal  to  the  corresponding 
component  of  attraction  at  0  reckoned  for  a  vanishing  circular  cavity 
with  0  as  centre  ;  the  limit  is  the  same  from  whichever  side  of  the 
surface  P  moves  up  to  0. 

33.  In  the  case  of  the  normal  attraction  component  Z,  the  subject 
of  integration  is  a-  (z-  £)  r~3,  and  the  corresponding  component  at  0 
is  represented  by  a  convergent  integral.  Thus 

D  =  Urn  Lim      <r  *  -  £)  r~*  dS. 


f  <r  (*  -  £)  r 

J 


The  most  important  terms  of  the  error  introduced  into  the  integral 
of  this  formula  by  putting  <TO  for  <r,  r  for  r,  and  dS0  for  dS,  are  of  the 
same  order  of  magnitude  as 


or 


f  *„(*- 


for  both  of  these  X  =  l,  /x  =  2,  and  therefore,  by  §31,  the  error  has 
a  zero  limit.     Accordingly  D  is  the  limit  of 


and  in  the  former  of  these  we  notice  that  z  is  of  the  order  of  smallness 
of  s2,  and  therefore  the  integral  of  the  same  order  as  I  a0s2r'~3dS0,  which 
has  X  =  0,  ft  =  1,  and  therefore  the  limit  zero.     Thus 
D  —  —  Lim  Lim  o-0  I  £  r'~3dS0, 

e-*-0<c-*0        J 

the  integral  being  now  taken  over  the  plane  circular  area  bounded  by 
r'  =  c,  3  =  0.  Now  if  du  represent  the  solid  angle  subtended  by 
dS0  at  P,  \£[r/~ad89  =  dv,  and  the  integral  is  ±fdto,  the  sign  being 
positive  if  £  is  positive,  negative  if  £  is  negative.  If  we  make  K-*-0 
before  e  -*-  0,  clearly  the  limit  of  fda>  is  27r,  and  so  we  get 


the  upper  sign  corresponding  to  £  positive.  So  the  limit  of  Z  differs 
from  the  value  of  Zi  at  0  by  27ro-0,  the  excess  of  the  former  over  the 
latter  corresponding  to  an  attraction  27ro-0  towards  the  surface;  the 
difference  between  the  limits  of  Z  as  P  approaches  0  from  different 

*  If  we  make  the  same  assumption  with  regard  to  a  as  in  the  footnote  of  §  32 
the  index  of  s  will  be  m  in  this  integral. 
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sides  of  the  surface  is  47ro-0  ,  and  the  arithmetic  mean  of  these  limits 
is  the  value  of  Z  at  0. 

34.     The  potential  at  P  of  a  double  sheet,  or  normally  magnetised 
shell,  of  strength  //.'  at  the  point  Q,  is  given  by 


where  \j/  is  the  angle  between  QP  and  the  normal  at  Q  drawn  in  the 
sense  for  which  ^  is  reckoned  positive.  The  error  introduced  into  the 
subject  of  integration  by  taking  ^  at  points  near  0  to  mean  the  angle 
between  QP  and  the  axis  of  z,  and  so  replacing  cos  $  by  -(z-£)  r~l, 
corresponds  to  dropping  a  factor  which  differs  from  unity  by  a  quantity 
of  the  order  of  s2,  and  the  integral  of  this  error  taken  over  the  cap  is 
one  for  which,  in  the  notation  of  §31,  A  =  1,  ft=  1,  and  therefore  has 
a  zero  limit.  Hence  the  potential  of  a  double  sheet  has,  for  purpose 
of  finding  D,  the  same  form  as  the  integral  investigated  in  the 
preceding  Article  ;  and  the  limit  of  V  as  P  approaches  0  from  the 
positive  side  of  the  sheet  exceeds  by  47r/x0'  the  limit  as  P  approaches 
0  from  the  negative  side. 

35.  The  potential  of  a  surface-distribution  of  gravitating  matter, 
whose  surface  density  is  free  from  such  peculiarities  as  would  render 
invalid  the  properties  already  established,  has  in  a  certain  sense  a 
space  differential  coefficient  in  any  direction  at  any  point  0  of  the 
surface  ;  this  is  not  a  differential  coefficient  as  generally  defined,  since 
it  is  a  limit  which  has  different  values  according  as  the  consecutive 
point  P  approaches  0  from  one  side  of  the  surface  or  from  the  other. 
It  is  to  be  noticed  that  the  existence  of  a  differential  coefficient  cannot 
be  inferred  from  the  physical  property  that  force  equals  gradient  of 
potential,  since  0  is  a  point  not  in  free  space,  but  in  the  gravitating 
matter. 

The  theorem  is  that 


Lim{(F0-T 

OP-^O  OP-*0 

where  Fp  is  the  component  of  the  force  at  P  resolved  along  the 
tangent  to  the  path  by  which  P  approaches  0. 

To  prove  this  we  must  show  that  if  rj  be  any  arbitrary  small 
quantity  we  can  always  choose  a  point  K  on  the  curve  by  which  P 
approaches  0  such  that  for  all  positions  of  P  on  the  curve  between  K 
and  0 

V^-  V, 


OP 

where  FQ  represents  Lim  Fp. 
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Let  us  regard  r)  as  the  sum  of  three  arbitrary  parts  r/u  >/2,  and  rj3. 
We  take  a  point  J  on  the  curve  between  P  and  0,  and  notice  that 

V  —  V 

v  o      Y  P     ™ 

OP  °~      OP 

and,  remembering  that 


where  F  is  the  tangential  force  and  ds  an  element  of  the  .  curve,  we 
apply  the  first  theorem  of  mean  value  and  so  get 


where  Q  is  some  point  on  the  curve  between  /  and  P. 
Thus 

V    -V  V    -V 

V  O  ^  P  Tl      _        ^  O  ^  J 


Since  F  is  definite  at  all  points  between  0  and  P,  and  has  a 
definite  limit  for  a  point  tending  to  coincidence  with  0,  there  is  a 
definite  superior  limit  to  the  absolute  value  of  F  for  the  points  of  the 
curve  lying  between  0  and  any  definite  point  K  ;  we  call  this  superior 
limit  M. 

Now  we  choose  K  so  near  to  0  that,  for  all  points  P  between  K 
and  0,  |  Fp  -  FQ  <  rji  and  therefore  also  |  FQ  —  FQ  <  ^  ;  this  we 
can  do  because  Fp  has  the  limit  FQ. 

We  next  take  P  anywhere  on  OK,  and  P  having  been  chosen,  we 

can  choose  a  point  L2  so  that  for  all  points  J  between  L2  and  0, 

VQ-  Vj\  <  OP  .  rj2)  this  being  possible  because  Vj  has  the  limit  VQ. 

And  a  point  Ls  can  be  chosen  so  that  for  all  points  J  between  L3 

and  0 


where  |  e  |  <  rjs/M.     We  now  take  J  to  be  between  0  and  the  nearer  of 
the  points  L2)  L3. 
Thus 

V     -  V         Tf>  V    —  V 

' 


r  V  —  V 
o       J 
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where  we  notice  that  |  F^  \  <  M.  The  modulus  of  the  first  expression 
in  square  brackets  is  less  than  173,  that  of  the  second  is  less  than  ijlt 
that  of  the  third  is  less  than  %,  hence  the  modulus  of  the  sum  of  the 
three  expressions  is  less  than  ^  -f  rj2  +  rj3  or  r).  Thus  we  have  been  able 
to  choose  Kso  that  for  all  points  P  on  the  curve  between  0  and  K 
V  -  V 

V  O  V  P  7j 

OP 

which  establishes  the  theorem. 


IX.     Volume  Integrals  through  regions  that  extend 
to  infinity. 

36.  The  integrals  to  which  we  have  so  far  been  devoting  most 
attention  are  those  whose  peculiarity  consists  in  the  subject  of 
integration  becoming  infinite  at  a  point  in  the  range.  Another  kind 
of  integral  requiring  special  study  occurs  frequently  in  mathematical 
physics,  namely,  a  volume  integral  taken  through  a  region  which 
extends  to  infinity. 

By  the  integral  ffdr  taken  through  all  space  outside  certain  finite 
closed  surfaces  Slt  S3)  etc.  is  meant  the  limit  of  the  integral  taken 
through  a  region  bounded  internally  by  $1,  $2,  etc.,  and  externally 
by  a  surface  B,  as  the  linear  dimensions  of  B  and  the  distances  of  all 
its  points  from  the  inner  boundaries  become  indefinitely  great,  provided 
such  a  limit  exists  and  is  independent  of  the  shape  of  B.  When  the 
limit  exists  and  is  independent  of  the  shape  of  B  the  integral  is  said 
to  be  convergent ;  if  the  limit  has  a  finite  value  which  is  not  inde- 
pendent of  the  shape  of  B,  the  integral  is  said  to  be  semi-convergent. 

The  following  is  the  chief  test  of  convergence.  If  we  measure  r 
from  some  fixed  origin,  and  if  f  is  such  that,  for  all  values  of  r  greater 
than  a  definite  length  a,  f  is  less  in  absolute  value  than  Mr~*  where  M 
is  a  constant  and  /*  >  3,  the  integral  is  convergent.  We  shall  prove 
this  by  shewing  that,  corresponding  to  any  arbitrarily  chosen  small 
quantity  o-,  there  can  always  be  found  a  closed  surface  0  surrounding 
0  and  all  the  surfaces  Si,  $2,  etc.,  such  that  all  closed  surfaces  t 
surrounding  6  have  the  property  that 


I/, 


fdr 


Take  for  the  surface  0  a  sphere  r--=iy  large  enough  to  surround  the 
sphere  r  -  a  and  all  the  inner  boundaries  of  the  region  ;  and  let  <o  be 
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the  distance  from  0  to  the  furthest  point  of  the  outer  boundary  t. 

Then 

rt 
fdr 


f  \f\dr,    <T  \f\dr, 

Je  Jo 

the  upper  limit  in  the  last  integral  being  the  sphere  r  =  o>.     Thus 


fdr 


<M 


/-o>  rot 

r'^dr,  <±irM  \   r*-*dr, 

J6  J-n 


,  a  positive  quantity, 


Hence  by  choosing  rj  greater  than  {47rM/(/x  -  3)  o-}^-3  we  get  a  surface  6 
satisfying  the  specified  condition ;  the  integral  is  accordingly  con- 
vergent. It  will  be  noticed  that  there  is  no  restriction  on  the  shape 
of  the  outer  boundary  t. 

Generally  speaking,  if  f  is  zero  at  infinity  of  an  order  higher  than 
r~*  the  integral  is  convergent ;  if  the  zero  is  just  of  the  order  r~3,  the 
integral  may  be  semi-convergent  or  divergent. 

37.  When  Green's  theorem  and  allied  theorems  are  applied  to 
volume  integrals   of  this  type,  the   outer   boundary  which  tends  to 
become  infinitely  large  must  not  be  left  out  of  account,  and  so  we 
have  limits  of  surface  integrals  which  are  spoken  of  as  integrals  over 
the  surface  infinity.     If  the  subject  ^  of  integration  is,  for  values  of 
r  greater  than  a  finite  length  a,  less  in  absolute  value  than  Mr~*9 
where  M  is  finite  and  /*>  2,  then  \tydS\  taken  over  the  sphere  r  =  w 
is  less  than  MvP~*Jfsxn  6  dO  d<f>,  which  has  the  limit  zero  for  o>  -*  oo . 
Thus  the  surface  integral  vanishes  if  ^  is  zero  at  infinity  of  a  higher 
order  than  r~z.     If  \j/  is  zero  of  the  order  r~2  the  limit  of  the  surface 
integral  may  be  different  for  different  shapes  of  B  ;  if  this  is  the  case 
there  is  of  course  corresponding  semi-convergence  of  one  of  the  volume 
integrals,  and  special  investigation  is  required. 

38.  The  differentiation  with  respect  to  a  parameter  of  a  volume 
integral  through  a  region  extending  to  infinity,  involving  as  it  does 
two  distinct  passages  to  limits,  requires  special  consideration.     Let 
us  consider  the  case  in  which  the  parameter  £  affects  the  subject  of 
integration,  but  does  not  affect  the  specification  of  the  inner  boundaries 
Siy  $2,  etc.     Let  us  suppose  that  df/dg  (or/')  exists  and  is  uniformly 
continuous  through  all  finite  portions  of  the  region  of  integration  for 
all  values  of  $  considered,  and  that  the  integral  of  /  is  convergent ; 


46 
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and  further  that,  for  all  values  of  £  considered  and  for  all  values  of  r 
greater  than  a,  there  is  an  inequality  \f  <Mr~^t  where  ^>3,  and 
Mt  /«,  and  a  are  constants  whose  values  do  not  depend  on  the  value 
of  £*.  Take  the  outer  boundary  to  be  the  sphere  r  =  w  ;  then 


•  •  .(28), 


=  Lim    Lim  M-  f  "  [f((  +  A  i  )  -/<*)}  *  -  /"/' 

A£-^0     <o-**>  L  A£J  / 

and  this,  by  the  theorem  of  mean  value, 

/M 
€C?T, 
—  ,    -    -      -    - 

where  «  =/'  (£  +  0  A  £)  —  f  (£)  and  1  >  0  >  0,  and  the  notation  implies 
that  first  w  -^  GO  and  afterwards  A  £  -*•  0.  If  we  can  shew  that  the 
subject  of  this  double  limit  can,  by  first  making  <o  -^-  GO  ,  and  afterwards 
taking  A  £  sufficiently  small,  be  made  less  than  any  assigned  small 
quantity  a-,  clearly  the  double  limit  will  be  zero. 

Now  since  /'  satisfies  the  conditions  of  the  theorem  of  §  36,  and 
moreover  in  such  a  way  that  M,  p.,  and  a  are  independent  of  £,  it  is 
clear  by  the  reasoning  of  that  Article  that,  for  all  values  of  £  considered 
and  therefore  in  particular  for  all  possible  values  of  £  +  6  A  £,  we  can 
choose  a  definite  length  77  such  that  for  all  values  of  w  greater  than  rj 

and 
are  both  less  than  any  assigned  small  quantity,  which  we  shall  take 


to  be  Jo-. 


and 


Hence 


Lim 


Liml  f 

w-*-oo    I  J  n 


JCT. 


And  rj  being  chosen  and  therefore  finite,  however  large,  the  uniform 
continuity  of/'  ensures  our  being  able  to  choose  a  value  of  A£  such 
that  for  it  and  for  all  smaller  values  j  e  is  less  than  an  arbitrary  small 
quantity  ;  this  small  quantity  we  choose  to  be  Jo-  T~l,  where  T  is  the 

r*j 


finite  volume      dr.     This  makes 


-Lim 


<  Jo-.     Thus 


/ 


<  (T. 


*  We  can  get  greater  generality  by  simply  requiring  that  the  integral  of/'  shall 
be  uniformly  convergent  for  the  contemplated  range  of  values  of  £ ;  but  it  seems 
better  not  to  introduce  into  the  text  the  idea  of  uniform  convergence,  especially  as 
there  are  additional  difficulties  in  the  proof. 
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Hence  the  double  limit  on  the  right-hand  side  of  (28)  is  zero,  and 
therefore  the  differentiation  of  the  integral  of  /  is  effected  by  the  rule 
of  differentiating  under  the  sign  of  integration. 

Differentiation  with  respect  to  a  parameter  which  affects  only  the 
specification  of  one  of  the  inner  boundaries,  say  Si,  clearly  gives  rise 
merely  to  a  surface  integral  over  Slm 

39.  Volume  integrals  through  regions  extending  to  infinity  occur 
in  electrical  theory  as  expressions  for  electrostatic  and  for  electro- 
dynamic  energy,  and  in  other  ways.  They  occur  in  the  theory  of 
gravitation  and  electrostatic  potential  in  proofs  of  the  important 
"theorems  of  uniqueness."  They  occur  in  Hydrodynamics  as 
representing  kinetic  energy,  and  "  impulse."  Differentiation  of  such 
integrals  is  employed  in  the  dynamical  theory  of  solid  bodies  moving 
through  an  infinitely  extended  liquid.  In  every  such  application  of 
these  integrals  it  is  necessary  to  make  sure  that  there  is  such 
convergence  as  will  render  the  formulae  valid. 
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